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ABSTRACT 


A  thermodynamic  model  has  been  developed  for  the  solubility  of  acid  gases  in 
alkanolamine  solutions.  The  model  is  based  on  the  extended  Debye-Huckel  theory  of 
electrolyte  solutions.  The  interaction  parameters  are  obtained  from  the  equilibrium  solubility 
data  for  the  pure  acid  gases  and  their  mixtures.  The  model  was  used  to  correlate  the  solubility 
of  carbon  dioxide,  hydrogen  sulfide,  and  their  mixtures  in  monoethanolamine,  diethanolamine 
and  diisopropanolamine  solutions,  in  the  temperature  range  of  25  to  120°C.  The  predicted 
partial  pressures  are  in  good  agreement  with  the  experimental  data.  However,  the  agreement 
for  the  partial  pressures  of  carbon  dioxide  over  3.5  N  DEA  solutions  is  not  satisfactory. 
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NOMENCLATURE 


a  -  Activity 

a  -  Ratio  of  activities  (amine  to  water) 
b  -  A  constant  in  Debye-Huckel  theory 
c  -  Concentration,  kmol/m3 
c  -  Total  number  of  moles,  mol 
d  -  Differential 
m  -  Molality,  mol/kg-solvent 
n  -  Number  of  moles,  mol 
u  -  Mole  ratio 

x  -  Mole  fraction  in  the  liquid  phase 

y  -  Mole  fraction  in  the  vapor  phase 

z  -  Electric  charge,  in  units  of  electronic  charge 

B  -  Array  of  interaction  parameters 

C  -  Deviations  in  predictions 

C  -  Number  of  components 

E  -  Errors  in  prediction 

E  -  Number  of  components 

G  -  Gibbs  free  energy,  J 

H  -  Henry’s  constant 

K  -  Equilibrium  constant 

P  -  Pressure,  Pa 

P  -  Array  of  predictions 

Q  -  Heat,  J 

R  -  Gas  constant,  J/mol  K 
S  -  Entropy,  J/K 
T  -  Temperature,  K 
U  -  Internal  energy,  J 

V  -  Volume,  m3 
W  -  Work,  J 

Y  -  Chemical  formula 


Greek  letters 

a  -  Stoichiometric  coefficient 

a  -  Mole  ratio  in  liquid 

(3  -  Binary  interaction  parameter 

7  -  Activity  coefficient 

A  -  Lagrange  multiplier 

q  -  Chemical  potential,  J/mol 

</>  -  Fugacity  coefficient 


IX 


A  -  A  small  perturbation 
II  -  Product  sign 
X  -  Summation  sign 

Subscripts 

bh  -  Acidic  dissociation  of  amine 
i  -  Component  i 
i  -  Dissociation  of  amine 
j  -  Element  j 
r  -  Reaction  r 
A  -  Carbamate  formation 
M  -  Carbamate  formation 
W  -  Water 
0  -  Solvent 

1C  -  Ionization  of  hydrogen  sulfide 
2C  -  Dissociation  of  bisulfide  ion 
1Y  -  Ionization  of  carbon  dioxide 
2Y  -  Dissociation  of  bicarbonate  ion 
+  -  Positively  charged  species 
-  -  Negatively  charged  species 

Superscripts 

0  -  Reference  state 
L  -  Liquid  phase 
V  -  Vapor  phase 

Symbol 

[...]-  Concentration,  kmol/m’ 
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I.  INTRODUCTION 


Natural  and  manufactured  gases  frequently  contain  objectionable  impurities  like 
hydrogen  sulfide  and  carbon  dioxide,  with  small  quantities  of  mercaptans,  carbon  monoxide, 
hydrocyanic  acid,  carbon  disulfide  and  carbonyl  sulfide.  Hydrogen  sulfide  is  particularly 
undesirable  on  account  of  its  toxicity  and  corrosiveness.  Removal  of  carbon  dioxide  is 
necessary  if  the  gas  is  to  be  liquefied  or  used  in  a  chemical  process. 

Aqueous  solutions  of  alkanolamines  are  the  most  widely  employed  solvents  for  the 
absorption  of  carbon  dioxide  and  hydrogen  sulfide.  The  use  of  alkanolamines  for  the  removal 
of  acid  gases  was  first  suggested  by  Bottoms  (1931)  who  discovered  that 

“certain  amines  and  amino  compounds  have  the  property  of  absorbing  carbon  dioxide 
and  hydrogen  sulfide,  forming  pyrolytically  unstable  carbonates  and  hydrosulfides 
which,  when  heated  to  a  temperature  above  50°C,  decompose  and  liberate  the  carbon 
dioxide  and  hydrogen  sulfide  and  release  the  free  base  in  its  original  form.” 

The  alkanolamines  are  characterized  by  the  presence  of  the  hydroxyl  group  -OH,  and 
the  amino  group  -NH2.  The  hydroxyl  group  reduces  the  vapor  pressure  of  the  amine  and 
increases  its  water  solubility,  while  the  amino  group  provides  the  necessary  reactivity  with  the 
acid  gases  (Kohl  and  Riesenfeld,  1974). 

Triethanolamine  (TEA)  was  the  first  to  be  used,  but  it  has  been  superseded  by 
monoethanolamine  (MEA)  and  diethanolamine  (DEA).  These  amines  have  a  higher  capacity, 
higher  reactivity  and  better  thermal  stability  than  TEA.  DEA  is  the  preferred  solvent  for  gases 
containing  carbonyl  sulfide  because  MEA  reacts  irreversibly  with  carbonyl  sulfide. 
Diisopropanolamine  (DIPA)  is  used  in  the  Adip  process  and  the  Sulfinol  process. 
Methyldiethanolamine  (MDEA)  is  a  selective  absorbent  for  hydrogen  sulfide  in  the  presence  of 
carbon  dioxide.  Diglycolamine  (DGA),  which  is  a  hydroxyaminoether  (and  not  a  true 
alkanolamine),  has  a  lower  vapor  pressure  than  MEA.  It  can  be  used  in  high  concentrations, 
resulting  in  higher  acid  gas  pickup,  lower  steam  consumption  and  a  decrease  in  the  amount  of 
solution  circulated.  Table  1  shows  the  structural  formulas,  molecular  weights  and  the  chemical 
names  of  these  amines. 

The  basic  flow  scheme  for  the  acid  gas  removal  process  is  shown  in  Figure  1.  The 
sour  gas  is  passed  upward  through  the  absorber,  countercurrent  to  the  amine  solution.  The 
sweetened  gas  leaves  the  top  of  the  absorber  and  is  dehydrated  before  being  put  to  use.  The 
rich  solution  from  the  bottom  of  the  absorber  is  heated  in  the  amine-amine  heat  exchanger 
and  fed  to  the  top  tray  of  the  stripping  column.  The  rich  solution  flows  down  the  stripping 
column  in  countercurrent  contact  with  the  vapor  generated  in  the  reboiler  and  is  stripped  of 
the  acid  gases.  Some  of  the  lean  solution  leaving  the  stripping  column  is  sent  to  the  reboiler. 
The  remaining  solution,  together  with  the  condensate  from  the  reboiler  is  passed  through  the 
amine-amine  heat  exchanger  and  a  cooler  before  being  fed  to  the  top  tray  of  the  absorber. 

The  acid  gas  leaving  the  stripping  column  is  cooled  to  condense  the  water  vapor.  The 
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condensate  is  separated  in  the  separator  and  is  fed  back  at  the  top  of  the  stripping  column. 

This  basic  scheme  may  be  modified  by  including  (1)  water  wash  of  the  purified  gas 
for  amine  recovery,  (2)  flashing  of  the  rich  solution  to  remove  dissolved  hydrocarbons  or 
(3)  split-stream  flow  to  reduce  the  steam  requirement. 

Over  the  past  fifty  years,  several  investigators  have  measured  the  equilibrium 
solubility  of  acid  gases  in  alkanolamine  solutions.  Kohl  and  Riesenfeld  (1974)  and  Lee  et  al. 
(1973)  have  presented  literature  reviews. 

The  rates  of  absorption  and  stripping  depend  upon  the  departure  from  equilibrium  in 
the  absorber  and  the  stripping  column.  Accurate  equilibrium  data  are  essential  for  the  proper 
design  of  the  gas  removal  process.  The  solubility  data  have  customarily  been  published  as  the 
plots  of  the  acid  gas  partial  pressure  against  the  mole  ratio  in  liquid  (acid  gas/amine),  at  a 
given  temperature  and  amine  concentration.  The  solubility  of  the  mixtures  is  also  expressed  in 
the  same  fashion,  with  the  mole  ratio  of  the  other  acid  gas  as  a  parameter. 

These  data  are  useful  in  many  ways.  They  can  be  used  to  calculate  the  acid  gas 
loadings  at  given  partial  pressures,  the  minimum  circulation  rate  of  the  amine  solution  and  the 
number  of  trays  in  the  absorber  and  the  stripping  column.  In  existing  plants,  these  data  may 
be  used  to  determine  how  close  the  absorber  is  to  equilibrium  at  the  lean  end  and  if  a  pinch 
condition  exists  at  the  rich  end.  As  the  amine  solution  flows  down  the  absorber,  its 
temperature  increases  because  of  the  heat  released  by  the  absorption  of  acid  gases.  Therefore, 
the  solubility  data  are  needed  for  many  different  temperatures.  The  effort  required  to  collect 
these  data  and  the  wide  range  of  process  conditions  make  a  computer-oriented  correlation  of 
the  data  highly  desirable.  This  correlation  would  be  useful  for  interpolating  the  solubility  data 
to  temperatures,  loadings  or  amine  concentrations  where  no  data  are  available. 

Simple  curve  fitting  is  unsuitable  because  it  provides  no  insight  into  the  physical 
basis  for  solubility  and  it  can  not  be  used  with  confidence  at  extremes  of  loadings  or 
pressures.  Such  a  model  would  be  inflexible  and  difficult  to  maintain  as  new  data  could  not 
be  easily  incorporated  into  it.  The  model  would  contain  a  large  number  of  meaningless 
coefficients.  The  nature  of  interpolation  is  also  not  obvious.  For  example,  it  is  not  obvious 
whether  the  same  functions  should  be  used  for  interpolating  over  temperatures  as  those  used 
for  interpolating  over  normalities. 

A  thermodynamic  model  based  on  the  physical  processes  occurring  in  the  system 
would  certainly  be  preferable  to  the  curve  fitting  approach.  This  model  would  be  able  to 
correlate  the  solubility  data  by  making  realistic  approximations.  Recently,  two  models  have 
been  proposed  to  correlate  the  equilibrium  data  for  acid  gas  /  alkanolamine  systems.  Before 
analyzing  these  models  in  detail,  it  is  necessary  to  understand  the  thermodynamic  criteria  for 
equilibrium  in  multicomponent,  multiphase  systems. 
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Table  1.  Alkanolamines 


(1)  Monoethanolamine 
Abbreviation  :  MEA 
Chemical  formula  :  H2NCH2CH2OH 
IUPAC  name  :  2-aminoethanol 
Molecular  weight  61.09 

(2)  Diethanolamine 
Abbreviation  :  DEA 

Chemical  formula  :  HN(CH2CH2OH)2 
IUPAC  name  :  2,2'-dihydroxydiethylamine 
Molecular  weight  105.14 

(3)  Triethanolamine 
Abbreviation  :  TEA 

Chemical  formula  :  (HOCH2CH2)3N 
IUPAC  name  :  2,2',2"-trihydroxytriethylamine 
Molecular  weight  149.19 

(4)  Diisopropanolamine 
Abbreviation  :  DIPA 

Chemical  formula  :  [CH3CH(OH)CH2]2NH 
IUPAC  name  :  2,2'-dihydroxydipropylamine 
Molecular  weight  133.19 

(5)  Diglycolamine 
Abbreviation  :  DGA 

Chemical  formula  :  HOCH2CH2OCH2CH2NH2 
IUPAC  name  :  2,2'-hydroxyaminoethylether 
Molecular  weight  105.14 

(6)  Methyldiethanolamine 
Abbreviation  :  MDEA 
Chemical  formula  :  (HOCH2CH2)2NCH3 
IUPAC  name  :  2,2'-dihydroxydiethylmethylamine 
Molecular  weight  119.17 


. 


4 


Acid  gas 


Figure  1.  Acid  Gas  Removal  Process 


II.  THERMODYNAMIC  FRAMEWORK 


The  criteria  for  equilibrium  can  be  derived  from  the  laws  of  thermodynamics.  These 
criteria  allow  us  to  examine  the  assumptions  implicit  in  various  solubility  models  and  point 
the  way  to  areas  where  improvements  can  be  made. 

The  second  law  of  thermodynamics  states  that 

dS  +  dSsurr  ^  0 

If  dQ  is  the  heat  added  to  the  system, 
dSsurr  =  -  dQ  /  T 

The  first  law  of  thermodynamics  for  closed  systems  is 
dQ  =  dU  -  dW 

If  only  mechanical  work  is  performed  on  the  system, 

dW  =  -  P  dV  —  (2.1) 

These  equations  can  be  combined  to  give 
dS  -  (dU  +  P  dV)  /  T  ^  0 
This  can  be  stated  in  the  form 
dU  +  P  dV  -  T  dS  ^  0 

as  T  is  never  less  than  zero.  The  Gibbs  free  energy,  G,  is  defined  by 
G  =  U  +  PV-TS 
At  a  constant  temperature  and  pressure, 
dG  =  dU  +  P  dV  -  T  dS 

The  second  law  can  now  be  stated  in  terms  of  the  Gibbs  free  energy, 

(dG)  ^  0  at  constant  T  and  P  —  (2.2) 

Any  change  occurring  within  the  system  has  to  satisfy  this  condition.  A  state  of  equilibrium  is 
attained  when  G  reaches  its  minimum  value.  Further  change  is  impossible  as  this  change 
would  not  satisfy  Equation  (2.2). 

For  a  p-phase  system  consisting  of  C  components,  the  fundamental  equation  can  be 
written  in  the  form 

G  =  G  (  T,  P,  nT  )  i  =  1,  .  .  .  ,  C;  n  =  1,  .  .  .  ,  p 

where  nT  is  the  number  of  moles  of  the  component  i  in  the  phase  tt.  The  chemical  potential 
IMn  is  defined  by 

MT  =  dG  /  9nT 

In  this  expression,  the  partial  derivative  is  taken  with  respect  to  nT,  keeping  the  temperature, 
pressure  and  the  number  of  moles  of  each  component  in  every  phase  (except  nT)  constant. 

The  Gibbs  free  energy  is  an  extensive  property  of  the  system. 

G  (  T,  P,  k  nT  )  =  k  G  (  T,  P,  nT  ) 

Applying  Euler’s  theorem,  we  have 
G  =  X  nT  mn 
Therefore, 
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dG  =  2  n,"  d/i"  +  2  mi"  dni" 

However,  differentiating  the  fundamental  equation, 
dG  =  X  Mi^  dn," 

These  two  equations  yield  the  Gibbs-Duhem  equation 
X  n,"  dMi"  =  0 

In  the  state  of  equilibrium,  G  is  minimum  and  the  variations  in  G,  caused  by  small  variations 
in  n,"  must  be  zero.  The  variations  dn,"  are  not  independent  because  they  have  to  satisfy 
additional  constraints  such  as  mass  balance,  charge  balance,  non-negativity,  etc..  The  solution 
to  this  constrained  minimization  problem  requires  Lagrange  multipliers. 

Let  u,j  be  the  moles  of  element  j  in  one  mole  of  component  i  and  let  q  be  the  total 
number  of  moles  of  element  j  present  in  the  system.  If  E  is  the  number  of  different  chemical 
elements  present  in  the  system,  the  mass  balance  can  be  written  as 

2  n,"  Uij  -  q  =  0  —  (2.3) 

The  charge  balance  can  be  expressed  in  this  form,  with  Ce+i  =  0  and  Ui.e+i  =  the  charge  on  the 
i-th  species. 

An  augmented  function  GA  can  be  defined  as 
Ga  =  G  +  2  Aj  (  2  n,"  Ujj  q  ) 

The  condition  of  equilibrium  now  becomes 
<5Ga  =  0 
Therefore, 

2  Mi77  dn,"  +  2  dAj  (  2  n,"  u.  -  c)  *  2  A,  (  2  Uy  dnj"  )  =  0 
And  thence, 

2  (  M'77  +  2  Aj  Uij  )  dns"  +  2  (  2  n;"  u,j  -  q  )  dAj  =  0 
As  dni"  and  dA,  are  independent  variables,  each  coefficient  in  the  above  equation  must  be 
zero. 

Mi77  +  2  Aj  Uij  =  0  —  (2.4) 

(and  2  n,"  Uy  -  q  =  0,  which  gives  back  the  constraint  equations).  These  equations  can  be 
solved  for  n,"  and  Aj  to  obtain  the  equilibrium  composition  and  the  Lagrange  multipliers. 

Traditionally,  thermodynamicists  have  avoided  the  use  of  artificial  variables,  A;,  by 
postulating  chemical  reactions.  The  reactions  are  expressed  as 

2  Oi"r  Yi"  =  0 

where  Yj"  is  the  chemical  formula  for  the  i-th  component.  The  superscript  7r  indicates  the 
phase  in  which  the  reaction  is  taking  place.  The  coefficients  a,Kr  are  called  the  stoichiometric 
coefficients  for  the  reaction  r.  By  convention,  atnr  are  positive  for  the  products  of  the  reaction 
and  negative  for  the  reactants.  These  postulated  reactions  allow  us  to  eliminate  Aj  from  the 
criterion  of  equilibrium. 

Any  postulated  chemical  equations  have  to  be  consistent  with  the  constraints 
mentioned  earlier.  The  coefficients  a ,"r  have  to  satisfy  the  equation 

2  uy  ai"r  =  0  --  (2.5) 

i.e.  every  postulated  chemical  equation  has  to  conserve  the  mass  of  each  element  and  the 
electric  charge.  Multiplying  Equation  (2.4)  by  a  ",  and  summing, 

2  a,n,  n,n  +  2  A)  (  2  Uy  a,\  )  =  o 
From  Equation  (2.5),  the  second  term  is  zero. 
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2  a" r  n"  =  0  —  (2.6) 

This  equation  expresses  the  condition  for  phase  and  reaction  equilibria  in  a  closed  system  at  a 
constant  temperature  and  pressure. 

Since  it  is  difficult  to  visualize  chemical  potentials,  the  criteria  for  equilibrium  are 
often  expressed  in  terms  of  activities,  defined  by 

H i"  =  +  RT  In  a, 77  —  (2.7) 

where  a, 77  is  the  activity  and  p?”  the  reference  value  of  the  chemical  potential.  The  reference 
state  is  the  one  of  unit  activity. 

The  condition  of  equilibrium  in  terms  of  activities  is 
2  a"T  (  /i?”  +  RT  In  aj77  )  =  0 

2  In  a, 77  =  -  2  (  /i?”  /  RT  )  —  (2.8) 

The  right  hand  side  depends  only  on  the  reference  states  chosen  and  is  a  constant  at  a 
prescribed  temperature. 

Equation  (2.8)  is  the  starting  point  for  most  studies  of  phase  and  reaction  equilibria. 
The  vapor-liquid  equilibrium  (VLE)  is  given  by 
Y,  (liquid)  =  Y;  (vapor) 

The  condition  for  equilibrium  is 

In  a,v  -  In  a,L  =  (  p?L  -  M,ov  )  /  RT 
a,v  /  a,L  =  exp  (  (  n?L  -  Mi°V  )  /  RT  ) 

The  right  hand  side  is  called  the  Henry’s  constant,  Hi. 

a,v  =  H,  a,L  —  (2.9) 

The  Henry’s  constant  depends  only  upon  the  reference  states.  This  equation  can  be  further 
modified  by  introducing  two  coefficients,  0,  and  y\.  The  fugacity  coefficient,  0,,  is  defined  by 
a,v  =  0,  y,  P  /  (  0,  y,  P  )° 

and 

0i  *  1  as  P  *  0. 

If  the  reference  state  is  chosen  so  that 
0,°  =  1,  y,°  =  1  and  P°  =  1, 
we  get 

a,v  =  0,  y(  P  —  (2.10) 

We  thus  get  a  rather  simple  expression  for  the  activity  if  the  reference  state  is  that  of  the 
perfect  pure  gas  at  the  system  temperature  and  unit  pressure.  This  state  is  only  hypothetical 
because  0,  will  be,  in  general,  different  from  unity  at  finite  pressures.  Also  the  expression 
av  =  0,  yi  P 

is  meaningful  only  if  P  is  measured  in  the  same  system  of  units  in  which  P°  is  unity.  For 
example,  if  P°  is  1  atm,  P  has  to  be  expressed  in  atmospheres. 

The  fugacity  coefficient  has  been  defined  in  such  a  way  that  for  an  ideal  mixture  of 
perfect  gases,  0,  =  1.  There  is  no  perfect  liquid  law  and  hence  no  preferred  definition  for  the 
activity  coefficient,  y\.  The  activity  coefficient  is  defined  in  several  ways,  one  of  which  is 
ajL  =  71  m,  /  7i°  m,° 

And  when  2  m,  —  0, 

7i  —  1  for  solutes 


and 
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yo  —  1  for  the  solvent. 

Other  units  of  concentration  such  as  molarity,  normality  or  mole  fraction  can  be 
used  in  place  of  molality.  These  lead  to  activity  coefficients  different  from  those  defined  above. 
It  is  clear  from  the  definition  that  the  activity  coefficients  tend  to  unity  in  an  infinitely  dilute 
solution.  The  reference  state  is  that  of  a  hypothetical  infinitely  dilute  solution  of  unit  molality. 
This  leads  to  a  simple  expression 

a,L  =  yi  m,  —  (2.1 1) 

It  must  be  emphasized  that  mi  is  in  the  same  set  of  units  in  which  mj°  is  unity.  In  Equation 
(2.11),  mi  is  a  dimensionless  number. 

With  these  definitions,  the  VLE  relation  becomes 
0,  yi  P  =  Hi  yi  m,  —  (2.12) 

By  convention,  the  activity  coefficient  of  water  has  been  based  on  its  mole  fraction  rather  than 
molality.  The  reference  state  is  that  of  pure  water  at  the  system  temperature  and  pressure.  The 
VLE  relation  for  water  is 

<£w  yw  P  =  Hw  yw  xw  —  (2.13) 

The  mole  fraction  of  water  is  usually  close  to  unity  so  that  yw  =  1.  For  liquid  water  in 
equilibrium  with  its  vapor,  xw  =  1  and  yw  =  1.  At  low  pressures,  the  fugacity  coefficient  is 
also  close  to  unity.  Thus  the  Henry’s  constant  of  water  is  equal  to  its  vapor  pressure.  This 
conclusion  is  valid  if  the  temperature  is  less  than  150°C  and  if  the  mole  fraction  of  water  in 
the  liquid  phase  is  close  to  unity. 

For  aqueous  phase  chemical  reactions,  the  condition  for  equilibrium  is 

2  ofir  In  a,  =  -  2  air  (  ju,0  /  RT  ) 
or 

n  (  a,  )  **  air  =  exp  [  -  X  air  (  /  RT  )  ] 

The  right  hand  side  is  called  the  equilibrium  constant,  Kr  for  the  reaction  r.  Like  Hi,  this 
dimensionless  constant  depends  only  upon  the  reference  states. 

Kr  =  n  (  a,  )  **  air  —  (2.14) 

where 

a,  =  yi  mi  for  solutes 
and 

aw  =  Xw  for  water. 

Equations  (2.12),  (2.13)  and  (2.14)  describe  the  conditions  for  equilibrium  in  an 
aqueous  solution  in  contact  with  its  vapor.  To  be  able  to  use  these  equations,  we  need 
information  about  (/>,  and  y\  for  all  the  components  present  in  the  system,  as  well  as  the 
Henry’s  constants  and  the  equilibrium  constants.  When  this  information  is  available  from 
experiments,  it  can  be  used  to  predict  the  equilibrium  composition.  In  absence  of  detailed 
information,  we  have  to  construct  simple  models  which  will  help  us  extract  this  information 
from  the  incomplete  data. 


. 


III.  SOLUBILITY  MODELS 


The  first  attempt  to  correlate  the  solubility  data  for  the  CO2  /  amine  /  water  system 
was  made  by  Mason  and  Dodge  (1936).  They  found  equations 

“which  represented  the  data  fairly  well  over  a  limited  range  of  conditions  but  all 
attempts  to  broaden  the  applicability  resulted  in  such  complexity  that  it  defeated  the 
purpose  for  which  the  equation  was  intended.” 

Although  Mason  and  Dodge  do  not  describe  the  method  they  used,  it  is  evident  that  they 
adopted  a  curve  fitting  approach.  The  reactions  between  alkanolamines  and  carbon  dioxide 
had  not  been  studied  at  that  time. 

Van  Krevelen  et  al.  (1949)  developed  a  method  for  predicting  the  vapor  pressures  of 
hydrogen  sulfide  and  ammonia  over  aqueous  solutions.  Their  scheme  was  later  used  by 
Danckwerts  and  McNeil  (1967)  for  the  CO2  /  amine  /  water  system.  Kent  and  Eisenberg 
(1975)  used  a  modification  of  this  method  to  correlate  the  equilibrium  solubilities  in  CO2  / 

H2S  /  amine  /  water  system. 

Another  model  developed  from  the  work  of  Atwood  et  al.  (1957)  who  proposed  a 
method  for  calculating  equilibrium  composition  in  the  H2S  /  amine  /  water  system.  Their 
method  was  generalized  by  Klyamer  et  al.  (1973)  and  applied  to  the  CO2  /  H2S  /  amine  / 
water  system.  A  modification  of  this  method  was  suggested  by  Nasir  (1975).  Edwards  et  al. 
(1975)  used  the  extended  Debye-Huckel  theory  of  electrolyte  solutions  to  calculate  the 
liquid-phase  activity  coefficients.  They  also  devised  a  method  to  obtain  the  parameters  required 
in  their  model.  Recently,  Beutier  and  Renon  (1978)  have  improved  Edwards’  method  by  a 
more  accurate  representation  of  activity  coefficients. 

All  these  models  are  based  on  the  chemical  reaction  equilibria  in  the  liquid  phase. 
The  main  reactions  occurring  in  the  CO2  /  H2S  /  amine  /  water  system  are  as  follows: 


Ionization  of  water: 

FLO  =  H+  +  OfT  —(3.1) 

Dissociation  of  hydrogen  sulfide: 

H:S  =  Fl+  +  HS'  —  (3.2) 

Dissociation  of  carbon  dioxide: 

H20  +  CO:  =  H+  +  HCO3'  —  (3.3) 

Dissociation  of  alkanolamine: 

H.O  +  RR'NH  =  RR'NH2+  +  OH-  —  (3.4) 

Formation  of  carbamate: 

RR'NH  +  C02  =  RR'NCOO'  +  H+  —  (3.5) 

Dissociation  of  bisulfide  ion: 

HS'  =  H+  +  S=  —  (3.6) 

Dissociation  of  bicarbonate  ion: 

HCO3'  =  H+  +  C03=  —  (3.7) 
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In  these  equations  RR'NH  is  the  chemical  formula  of  the  alkanolamine.  Tertiary 
amines  (TEA  and  MDEA)  lack  the  extra  hydrogen  atom  and  do  not  form  carbamates  by 
reaction  5.  In  addition  to  the  reactions  1  to  7,  other  reactions  may  take  place  in  the  solution. 
The  amines  react  with  carbon  dioxide  to  form  heterocyclic  compounds.  The  amines  may  react 
with  impurities  in  the  sour  gas  such  as  carbonyl  sulfide  and  carbon  disulfide.  These  reactions 
have  little  effect  on  the  acid  gas  equilibria  as  only  a  small  fraction  of  the  reactants  is  used  up 
in  these  reactions.  The  alkanolamines  have  two  active  groups,  viz.  -NH2  and  -OH.  On  account 
of  this  bifunctional  nature  of  alkanolamines,  the  following  reactions  are  possible: 

Ri-OH  +  OH'  =  Ri-O  +  H20 
R,-0  +  C02  =  Ri-O-COO' 

These  reactions  are  excluded  from  consideration  because  the  amount  of  products  is  negligibly 
small  (Jensen  et  ai,  1954).  Astarita  et  al.  (1964)  report  that  these  reactions  take  place  only  in 
a  basic  solution  where  pH  is  11  or  greater;  these  reactions  can  be  neglected  as  the  pH  of  even 
a  slightly  carbonated  solution  is  less  than  10. 

The  chemical  reactions  1  to  7  are  accompanied  by  the  vapor-liquid  equilibria  of  the 
volatile  species. 

C02  (1)  =  C02  (g)  —  (3.8) 

H2S  (1)  =  H2S  (g)  -  (3.9) 

H20  (1)  =  H20  (g)  —  (3.10) 

Free  amine  may  also  be  present  in  the  vapor  phase.  Under  the  conditions  of  interest, 
the  vapor  pressures  of  the  amines  are  very  small  and  we  can  assume  that  the  amine  is  present 
only  in  the  liquid  phase.  The  thermodynamic  framework  developed  in  the  previous  chapter 
enables  us  to  write  down  the  equilibrium  relations  among  the  reacting  species.  These  equations 


are  as  follows: 

Kw  =  y(H*)  y(OH')  m(H+)  m(OH")  /  xw  --  (3.11) 

K,c  =  y(H+)  y(HS~)  mfH*)  m(HS~)  /  y(H2S)  m(H2S)  —  (3.12) 

K1y  =  y(H+)  y(HCOO  m(H+)  m(HC03l  /  y(C02)  m(C02)  xw  —  (3.13) 

K,  =  yIRR'NHj*)  y(OH^  m(RR'NH2+)  m(OH~)  /  y(RR'NH)  m(RR'NH)  xw  --  (3.14) 
Ka  =  y(H+)  y(RR'NCOO')  mfH*)  m(RR'NCOO)  /  y(RR'NH)  y(C02)  m(RR'NH) 
m(C02)  —(3.15) 

K2C  =  y(H+)  y(S=)  m(H+)  m(S=)  /  y(HS~)  m(HS~)  —  (3.16) 

K2Y  =  y(H+)  y(C03=)  m(COf)  /  y(HCO 3~)  m(HC03~)  —  (3.17) 

0( C02)  y(C02)  P  =  H(C02)  y(C02)  m(C02)  —  (3.18) 

0(H2S)  y(H2S)  P  =  H(H2S)  y(H2S)  m(H2S)  —  (3.19) 

(f> w  yw  P  :  (Pw)vp  xw  —  (3.20) 

The  reference  states  have  been  defined  in  Chapter  II.  Equations  (3.11)  to  (3.19) 
describe  the  rigorous  thermodynamic  criteria  of  equilibrium  in  the  CO2  /  H2S  /  amine  /  water 
system.  Equation  (3.20)  for  water  is  based  on  the  assumption  that  </>w  =  1  for  pure  water  at 
its  vapor  pressure. 

The  mass  balance  and  the  electric  charge  balance  are  given  by  the  following 


equations: 

mA  =  m(RR'NH)  +  m(RR'NH2+)  +  m(RR'NCOCT)  —  (3.21) 

mA  «(C02)  =  m(C02)  +  m(HC03_)  +  m(C03=)  +  m(RR'NCOO')  —  (3.22) 
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mA  a(H2S)  =  m(H2S)  +  m(HS^  +  m(S  )  —  (3.23) 

rrtfH*)  +  mCRR'NH^  =  m(OH")  +  mens')  +  m(HCOO  +  m(RR'NCOO')  +  2  m(S=)  + 

2  m(C03=)  —  (3.24) 

Here  <*(C02)  and  a(H2S)  are  the  mole  ratios  in  liquid  (carbon  to  nitrogen  and  sulfur  to 
nitrogen). 

Olofsson  and  Hepler  (1975)  give  a  set  of  best  values  for  the  thermodynamic 
equilibrium  constant,  Kw.  These  values  were  obtained  from  an  analysis  of  several  experimental 
investigations  and  are  consistent  with  other  thermodynamic  properties  such  as  the  standard 
enthalpy  of  ionization.  They  recommend  the  following  equation  for  Kw,  valid  upto  150°C. 

-  logx o  Kw  =  142  613.6/T  +  4229.195  logl0  T  -  9.7384  T  +  0.012  963  8  T2  -  1.150  68X10'5 

T3  +  4.602X1 0"9T*  -  8909.483  —  (3.25) 

Rao  and  Hepler  (1977)  have  similarly  analyzed  the  equilibrium  constants  and  the 
thermodynamics  of  ionization  of  aqueous  hydrogen  sulfide.  For  temperatures  between  0°C  and 
250°  C,  they  recommend  the  following  equation: 

log io  K,c  =  -  6045.2/T  +  106.67  -  37.744  logl0  T  —  (3.26) 

The  first  ionization  constant  of  carbonic  acid  was  measured  by  Read  (1975)  from 
25°C  to  250°C.  His  results  can  be  correlated  by  the  following  equation: 

In  Kiy  =  413.7285  +  4.431  786X10-2  T  -  67.3414  In  T  -  1.726  204X104/T  —  (3.27) 

The  most  reliable  measurements  of  the  dissociation  constants  of  the  alkanolamines 
are  those  of  Bates  and  coworkers.  Bates  and  Pinching  (1951)  recommend  the  following 
equation  for  the  acidic  dissociation  constant  of  MEA,  from  0°C  to  50°C. 

-  logx o  Kbh  =  2677.9 1/T  +  0.3869  +  0.000  427  7  T  —  (3.28) 

The  basic  dissociation  constant,  K,  can  be  obtained  from  Kbh  and  Kw  by  the  equation: 

K,  =  Kw/Kbh  —  (3.29) 

Bower  et  al.  (1962)  give  the  following  equation  for  the  acidic  dissociation  constant  of  DEA, 
from  0°C  to  50°C. 

-  logx o  Kbh  =  1830.15/T  +  4.0302  -  0.004  326  1  T  —  (3.30) 

Bates  and  Allen  (1960)  report  the  acidic  dissociation  constant  of  TEA,  from  0°C  to  50°C. 

-  logx o  Kbh  =  1341. 16/T  +  4.6252  -  0.0045666  T  —  (3.31) 

The  measurements  of  the  dissociation  constant  of  the  alkanolamines  have  been  reported  only 
for  the  temperature  range  of  0°C  to  50°C.  No  measurements  have  been  reported  for  DIPA, 
DGA  and  MDEA. 

Jensen  et  al.  (1954)  studied  the  reactions  between  carbon  dioxide  and  alkanolamines. 
They  report  the  equilibrium  constants  for  carbamate  formation  for  MEA  and  DEA  at  18°C. 
They  have  used  the  Bjerrum  expression, 

-  logx o  7  =  0.3  (Cion)0'3333 

to  calculate  the  equilibrium  constant.  McNeil  (1965)  studied  the  same  reaction,  but  used  a 
different  method  to  analyze  his  results.  Tertiary  amines  do  not  form  carbamates.  No  data  have 
been  reported  for  DIPA  and  DGA. 

For  moderate  degrees  of  carbonation,  the  carbamate  is  the  predominant  product  of 
the  reaction.  It  is  indeed  unfortunate  that  this  key  reaction  has  received  so  little  attention  from 
researchers.  At  present,  the  equilibrium  constants  for  carbamate  formation  are  available  only 
for  MEA  and  DEA  at  temperatures  of  18°C  and  20°C.  These  measurements  have  been 
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interpreted  on  the  basis  of  untested  models. 

The  remaining  two  reactions,  the  dissociation  of  the  bisulfide  and  bicarbonate  ions 
are  important  only  at  very  low  loadings  of  the  acid  gases.  The  second  ionization  constant  of 
hydrogen  sulfide  has  been  measured  by  Kryukov  et  al.  (1974).  They  give  the  following 
equation  for  the  dependence  of  K2C  on  temperature,  valid  from  25°C  to  250°C. 

-  log  10  K2C  =  2892.91 /T  +  4.448  -  0.002  514  T  —  (3.32) 

The  second  ionization  constant  for  the  carbonic  acid  has  been  measured  by  Harned  and 
Scholes  (1941),  Cuta  and  Strafelda  (1954)  and  Ryzhenko  (1963).  Their  measurements,  as 
reported  by  Kent  and  Eisenberg,  may  be  correlated  by  the  equation: 

In  K2Y  =  23.869  04  -  4.458  435X10-2  T  -  2.245  836  In  T  -  6433.628/T  —  (3.33) 

The  Henry’s  constant  for  carbon  dioxide  has  been  measured  by  Wiebe  and  Gaddy  (1940), 
Zel’venskii  (1937)  and  Ellis  (1959).  These  data  were  correlated  by  the  equation: 

In  H(C02)  =  -  4379.848/T  +  36.725  75  -  2.387  349X10"2  T  —  (3.34) 

Lee  and  Mather  (1978)  have  measured  the  Henry’s  constant  for  hydrogen  sulfide.  They  have 
also  reviewed  the  earlier  measurements  of  the  Henry’s  constant.  Their  data  were  correlated  by 
the  equation: 

In  H(H2S)  =  -  4788.603/T  +  38.4457  -  2.897  821X10'2  T  —  (3.35) 


1.  Kent  and  Eisenberg  Model 

Van  Krevelen  et  al.  used  the  following  method  for  calculating  the  vapor  pressures  of 
ammonia  and  hydrogen  sulfide  over  aqueous  solutions. 

Instead  of  the  thermodynamic  equilibrium  constant  and  the  Henry’s  constant,  van 
Krevelen  et  al.  used  the  pseudo-equilibrium  constants  defined  as  follows: 

K'  =  K  7(H2S)  7(NH3)  /  7(HS“)  7(NH4+)  =  m(HS~)  m(NH4+)  /  m(H2S)  m(NH3)  —  (3.36) 
H'(H:S)  -  H(H2S)  7(H2S)  /  4>(H2S)  =  y(H2S)  P  /  m(H2S)  —  (3.37) 

Van  Krevelen  et  al.  have  used  molarity  instead  of  molality  as  the  measure  of 
concentration.  The  use  of  molality  does  not,  however,  alter  the  essential  features  of  their 
method.  Van  Krevelen  et  al.  proposed  equations  for  relating  these  pseudo-equilibrium  constants 
to  the  ionic  strength  of  the  solution.  The  ionic  strength  is  defined  as: 

H  =  0.5  2  m,  zi2  —  (3.38) 

It  has  been  claimed  that  the  dependence  of  Ki'c  and  H'(H2S)  on  the  ionic  strength  is  linear. 

In  their  work,  van  Krevelen  et  al.  have  ignored  the  second  ionization  of  hydrogen  sulfide. 

They  have  also  neglected  the  concentrations  of  the  hydrogen  and  hydroxyl  ions. 

Danckwerts  and  McNeil  used  this  approach  to  predict  the  equilibrium  partial 
pressure  of  carbon  dioxide  over  carbonated  amine  solutions.  The  central  feature  of  this  model 
is  the  use  of  pseudo-equilibrium  constants  and  their  dependence  on  the  ionic  strength.  The 
ratio  of  the  pseudo-equilibrium  constant  at  a  certain  ionic  strength  to  that  at  zero  ionic 
strength  has  been  termed  the  “ionic  characterization  factor”.  The  most  serious  limitation  of 
this  model  is  that  the  ionic  strength  alone  is  not  sufficient  to  determine  the  ionic 
characterization  factors.  Secondly,  the  ionic  characterization  factors  may  not  be  linearly  related 
to  the  ionic  strength.  Thirdly,  the  dependence  of  these  factors  on  the  ionic  strength  can  not  be 
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determined  without  making  additional  simplifications.  For  example,  several  species  may  have 
to  be  neglected  from  the  balance  equations. 

Kent  and  Eisenberg  tried  to  correlate  the  solubility  data  using  the  published 
equilibrium  constants  without  using  the  ionic  characterization  factors.  They  found  that  the 
predicted  values  deviate  substantially  from  the  experimental  values.  Instead  of  using  the  ionic 
characterization  factors,  they  correlated  the  data  by  treating  K,  and  KA  as  variables.  They 
accepted  the  published  values  of  other  equilibrium  constants  and  found  K,  and  KA  by  fitting 
the  data  for  15.3  weight  percent  MEA  solution  and  20.5  weight  percent  DEA  solution.  The 
Kent  and  Eisenberg  model  is  based  on  a  set  of  reactions  that  is  chemically  equivalent  to 


reactions  (3.1)  to  (3.7).  The  model  equations  are  listed  below: 

Reactions: 

RR'NH2+  =  RR'NH  +  H+  —  (3.39) 

RR'NCOO”  +  H20  =  RR'NH  +  HC03~  —  (3.40) 

H20  +  C02  =  H+  +  HCCV  —  (3.41) 

H,0  =  H+  +  OH”  —  (3.42) 

HCCV  =  H4  +  C(V  —  (3.43) 

H2S  =  H4  +  HS"  —  (3.44) 

HS”  =  H4  +  S=  —  (3.45) 

C02  (1)  =  C02  (g)  -  (3.46) 

H2S  (1)  =  H2S  (g)  -  (3.47) 

Equilibrium  relations: 

K,  =  [H4]  [RR'NH]  /  [RR'NH,4]  —  (3.48) 

K2  =  [RR'NH]  [HC(V|  /  [RR'NCOO!  [H,0]  —  (3.49) 

K3  =  [H4]  [HCOjI  /  [C02]  [H20]  —  (3.50) 

K4  =  [H4]  [OH!  /  [H20]  —  (3.51) 

K5  =  [H!  [CO!  /  [HCO3!  —  (3.52) 

K6  =  [H!  [HS!  /  [H2S]  —  (3.53) 

K 7  =  [H!  [S!  /  [HS!  --  (3.54) 

H(C02)  =  P(C02)  /  [CO!  —  (3.55) 

H(H2S)  =  P(H2S)  /  [H2S]  —  (3.56) 

Balance  equations: 

[H!  +  [RR'NH,!  =  [OH!  +  [RR'NCOO!  +  [HCOj!  +  [HS!  +  2  [C03=]  +  2  [S=]  — 

(3.57) 

cA  =  [RR'NH]  +  [RR'NH,!  +  [RR'NCOO!  —  (3.58) 

cA  a(C02)  =  [CO,]  +  [HCO,!  +  [RR'NCOO!  +  [C03=]  —  (3.59) 

cA  a(H,S)  =  [H,S]  +  [HS!  +  [S=]  —  (3.60) 


In  these  equations,  the  quantities  in  the  square  brackets  are  the  molarities  of  various 
species,  with  the  exception  of  [H,0],  which  is  the  mole  fraction  of  water. 

It  would  be  incorrect  to  use  equations  like  (3.50)  with  K3  obtained  from  literature. 
The  activity  coefficients  even  in  moderately  charged  solutions  are  of  the  order  of  0.5.  The 
ratio  of  concentrations  in  such  solution  would  be  four  times  the  ratio  given  by  Equation 
(3.50).  It  is  clear  that  even  if  the  Kent  and  Eisenberg  model  succeeds  in  predicting  the  correct 
partial  pressures,  it  would  fail  to  predict  the  concentrations  of  the  ionic  species. 
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The  VLE  relations  are  similarly  inaccurate  as  they  do  not  include  the  fugacity 
coefficients.  If  the  acid  gas  partial  pressure  is  a  few  MPa,  the  fugacity  coefficient  may  be  far 
from  unity. 

To  improve  the  predictions  made  by  this  method,  Kent  and  Eisenberg  treated  Ki  and 
K;  as  adjustable  parameters  and  determined  them  by  fitting  the  model  to  experimental  data. 
They  were  able  to  obtain  excellent  fits  by  using  this  approach.  They  also  discovered  that  the 
constants  Ki  and  K2  obtained  in  this  manner  show  an  Arrhenius  dependence  on  temperature. 

Kent  and  Eisenberg  obtained  their  constants  by  fitting  the  data  for  2.5  N  MEA 
solutions.  Nasir  used  these  constants  to  predict  the  acid  gas  partial  pressures  over  5.0  N  MEA 
solutions.  He  found  that  the  errors  in  partial  pressures  are  often  in  excess  of  100  percent.  The 
Kent  and  Eisenberg  model  was  also  unable  to  predict  the  partial  pressures  accurately  when  the 
acid  gas  loadings  were  were  either  very  low  or  very  high. 

This  procedure  of  lumping  all  non-idealities  into  two  adjustable  parameters  is 
unsatisfactory  for  many  reasons.  The  thermodynamic  rigor  is  lost  if  the  experimentally 
measured  values  are  disregarded.  The  constants  Ki  and  K2  depend  on  the  normality  of  amine 
solution.  The  predictions  are  accurate  only  in  a  narrow  range  of  loadings  (0.2  to  0.7).  The 
constants  Ki  and  K2  are  obtained  by  fitting  the  pure  component  data;  i.e.  the  data  for  one 
acid  gas  and  amine  solution.  When  these  constants  are  used  to  predict  the  partial  pressures 
with  both  acid  gases  present,  the  predictions  are  highly  inaccurate.  The  attractive  feature  of 
this  model  is  that  the  model  equations  are  amenable  to  algebraic  manipulations,  thus 
permitting  hand-calculation  of  the  equilibria. 


2.  Klyamer  et  al  Model 

Atwood  et  al.  proposed  a  method  of  calculation  of  equilibria  in  the  H2S  /  amine  / 
water  system.  Their  scheme  was  used  by  Klyamer  and  Kolesnikova  (1972)  to  predict  the 
equilibria  in  the  CO2  /  amine  /  water  system.  Klyamer  et  al.  generalized  the  scheme  to  make 
it  applicable  to  the  CO2  /  H2S  /  amine  /  water  system.  The  central  feature  of  this  model  is 
the  use  of  “mean  ionic  activity  coefficient”.  The  activity  coefficients  of  all  ionic  species  are 
assumed  to  be  equal. 

This  assumption  has  some  basis  in  theory.  In  the  case  of  a  uni-univalent  electrolyte, 
the  assumption  is  certainly  a  valid  one.  Consider  the  reaction: 

MX  =  M+  +  X"  —  (3.61) 

Here  the  individual  activity  coefficients  7(M*)  and  y(X~)  (usually  written  as  7+  and  7_)  are 
not  measurable  properties.  By  virtue  of  the  charge  balance,  concentrations  of  the  anion  and 
the  cation  are  equal  to  each  other  and  can  not  be  varied  independently.  This  means  that  7+ 
and  7_  are  experimentally  inaccessible.  Only  the  product  7+7-  can  be  measured.  For  a 
uni-univalent  electrolyte,  the  mean  ionic  activity  coefficient  is  defined  by 

7  =  (7+7-)05  —  (3.62) 

The  quantity  7  is  measurable  and  is  sufficient  to  characterize  the  equilibria.  For  a  single 
electrolyte,  the  assumption  of  equal  activity  coefficients  is  tautologous.  For  a  mixture  of 
electrolytes,  the  Debye-Huckel  theory  gives  an  expression  for  the  activity  coefficients.  The 
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Debye-Huckel  theory  is  applicable  only  at  low  ionic  strengths.  The  electrolyte  solution  is 
assumed  to  be  a  collection  of  electric  charges  embedded  in  a  dielectric  medium.  The  resulting 
expression. 

In  7i  =  -  A  z,2  m0'5  —  (3.63) 

shows  that  the  ionic  coefficients  are  equal  in  a  dilute  solution  of  uni-univalent  electrolytes. 

However,  there  is  no  reason  to  believe  that  the  ionic  activity  coefficients  will  be 
equal  in  a  solution  of  high  ionic  strength.  The  assumption  is  still  good  if  only  one  anion  and 
one  cation  are  present  in  significant  amount.  This  is  generally  not  the  case  for  the  C02  /  H2S 
/  amine  /  water  system. 

The  Klyamer  et  al.  model  is  also  based  on  a  set  of  reactions  that  is  chemically 
equivalent  to  the  reactions  1  to  7.  These  reactions  are  not  the  same  as  those  used  by  Kent 
and  Eisenberg,  but  are  chemically  equivalent.  The  model  equations  are  presented  below. 


Reactions: 

RR'NH  +  H20  =  RR'NH2+  +  OH'  —  (3.64) 

H2S  =  H+  +  HS~  —  (3.65) 

H20  =  H+  +  OH'  --  (3.66) 

HS'  =  H+  +  S=  —  (3.67) 

2  RR'NH  +  CO:  =  RR'NH2+  +  RR'NCOO'  —  (3.68) 

H20  +  CO:  =  H+  +  HCO3'  —  (3.69) 

HC03'  =  H+  +  C03=  —  (3.70) 

Equilibrium  relations: 

K,  =  y2  m(RR'NH2+)  m(OH')  /  a  a2  m(RR'NH)  —  (3.71) 

Kw  =  y2  m(H")  m(OH~)  /  a  --  (3.72) 

Kic  =  y2  m(H+)  /  m(H2S)  —  (3.73) 

K2C  =  7  m(H+)  m(S=)  /  m(HS~)  —  (3.74) 

H'(H2S)  =  m(H2S)  /  p(H2S)  --  (3.75) 

Km  =  72  m(RR'NH2+)  m(RR'NCOO")  /  a2  a2  m2(RR'NH)  p(C02)  —  (3.76) 

Kiy  =  72  m(H+)  m(HC03l  /  a  m(C02)  —  (3.77) 

K2Y  =  7  m(H")  m(C03=)  /  m(HC03~)  —  (3.78) 

H'(C02)  =  m(C02)  /  p(C02)  —  (3.79) 


The  balance  equations  are  identical  to  Equations  (3.21)  to  (3.24).  In  these  equations,  7  is  the 
mean  ionic  activity  coefficient,  a  is  the  mole  fraction  of  water  and  (aa)  is  the  activity 
coefficient  of  the  free  amine.  H'(H2S)  and  H'(C02)  are  the  inverted  Henry’s  constants.  The 
mean  ionic  activity  coefficient  is  determined  from  a  plot  given  by  Atwood  et  al.  where  it  is 
related  to  the  ionic  strength  of  the  solution.  This  plot  was  obtained  by  fitting  the  data  for 
H2S  /  amine  /  water  system  and  smoothing.  Klyamer  et  al.  assume  that  the  relation  between 
the  mean  ionic  activity  coefficient  and  the  ionic  strength  is  the  same  in  the  CO:  /  H2S  / 
amine  /  water  system.  The  ratio  “a”  is  related  to  the  molality  of  of  free  alkanolamine.  This 
functional  relation  was  estimated  by  Atwood  et  al.  from  the  data  of  Sivertz  et  al.  (1940). 

The  Klyamer  et  al.  model  also  assumes  that  the  activity  coefficients  of  the  free  H2S 
and  free  C02  are  equal  to  unity.  The  fugacity  coefficients  are  also  assumed  to  be  equal  to 
unity. 

The  equality  of  activity  coefficients  makes  this  model  suitable  for  analytical  solution. 
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Simple  expressions  for  the  partial  pressures  may  be  derived  from  these  equations.  These 
expressions  have  been  derived  in  Appendix  A. 

The  Klyamer  et  al.  model,  with  7  =  1  and  dux  =  1,  is  algebraically  equivalent  to  the 
Kent  and  Eisenberg  model.  If  these  substitutions  are  made  in  the  expressions  for  the  Klyamer 
et  al.  model,  the  corresponding  expressions  for  the  Kent  and  Eisenberg  model  are  obtained. 

Nasir  made  two  improvements  in  the  Klyamer  et  al.  model.  He  incorporated  into  the 
model  the  fugacity  coefficients  based  on  the  Redlich-Kwong  equation  of  state.  He  made  7,  the 
mean  ionic  activity  coefficient  a  function  of  the  initial  amine  concentration  and  the 
temperature.  These  changes  resulted  in  some  improvement  in  accuracy,  especially  at  high 
pressures. 


3.  Proposed  Model 

The  main  problem  in  the  calculation  of  VLE  in  multicomponent  systems  is  the 
estimation  of  the  fugacity  coefficients  and  the  activity  coefficients.  Kent  and  Eisenberg  avoid 
this  problem  by  lumping  all  non-idealities  into  two  adjustable  parameters.  Klyamer  et  al. 
assumed  that  all  ionic  activity  coefficients  are  equal  and  that  these  can  be  estimated  from  an 
empirical  correlation.  The  failure  of  both  models  indicates  the  need  for  a  better  method  to 
estimate  the  activity  coefficients  and  especially  at  high  pressures,  the  fugacity  coefficients  have 
to  be  estimated  accurately.  The  fugacity  coefficients  are  related  to  the  volumetric  behavior  of 
the  vapor  phase. 

In  </>,  =  S  (  Z,  -  1  )  d  In  P  --  (3.80) 

Here  Z,  =  P7/RT  and  7  =  the  partial  molar  volume  of  the  i-th  component.  We  estimate  the 
fugacity  coefficients  using  the  Peng-Robinson  equation  of  state.  This  equation  has  proved  to 
be  superior  to  the  original  R-K  equation  and  the  Soave  modification  of  the  R-K  equation. 

The  Peng-Robinson  equation  is  of  the  form 

P  =  (  RT/(v-b)  )  -  a(T)/(  v(v+b)  +  b(v-b)  )  —  (3.81) 

The  fugacity  coefficient  of  component  i  can  be  calculated  from  the  following  equation: 

In  0i  =  (b,/b)(Z-l)  -  In  (Z-B)  (1/V/B)(A/B)X 

((2  2  y,  a,j/a)-b,/b)  In  ((Z+2.414  B)/(Z-0.414  B))  —  (3.82) 

where 

A  =  aP/R2T2,  B  =  bP/RT,  Z  =  Pv/RT,  a  =  2  y.y,^,  b  -  2  y.bj  and  ay  =  (l-6ijXataj)0-5  — 
(3.83) 

In  Equation  (3.83),  <5ij  is  an  empirically  determined  binary  interaction  coefficient.  This 
coefficient  is  independent  of  temperature  except  for  water-containing  systems. 

The  constants  aj  and  b,  are  related  to  the  critical  properties  and  the  acentric  factor 
of  the  component  j.  The  subscript  j  has  been  omitted  in  the  following  equations  for 
convenience. 

a(T)  =  a(Tc)a(Tr,  <x>) 
b(T)  =  b(Tc) 

a(Tc)  =  0.45724  R2  (Tc)2  /  Pc 
b(Tc)  =  0.07780  R  Tc  /  Pc 


< 

' 


17 


<*0  5(Tr,  a>)  =  1  +  k[1  -  (Tr)0  5] 

k  =  0.37464  +  1.54226  w  -  0.26992  a,2  —  (3.84) 

The  compressibility  factor  Z  is  determined  by  solving  Equation  (3.81).  The  equation 
yields  one  or  three  real  roots  depending  upon  the  number  of  phases  in  the  system.  In  the 
two-phase  region,  the  largest  root  is  the  compressibility  factor  of  vapor.  The  constants  required 
for  solving  Equation  (3.81)  are  easily  determined  from  Equations  (3.83)  and  (3.84).  The 
Peng-Robinson  equation  is  a  cubic  equation  and  it  can  be  solved  analytically.  The  predictions 
made  by  the  P-R  equation  are  particularly  good  in  the  critical  region.  This  is  advantageous  as 
the  critical  temperatures  of  CO2  and  H2S  fall  in  the  range  of  interest. 

The  Debye-Huckel  expression  may  be  used  for  the  prediction  of  activity  coefficients 
in  very  dilute  solutions.  For  more  concentrated  solutions,  Guggenheim  (1935)  gives  the 
extended  Debye-Huckel  expression 

In  7k  =  -  A  (zk)2  n°'5  /  (1  +  bk  /I05)  +  2  2  /3kj  mj  —  (3.85) 

The  first  term  expresses  the  Debye-Huckel  law  and  represents  the  electrostatic  forces. 
The  second  term  takes  into  account  the  short  range  van  der  Waals  forces.  The  quantity  A  is 
related  to  the  dielectric  constant  of  the  solvent.  Usually  bk  may  be  considered  a  constant. 
Values  of  1.0,  1.2  or  1.5  are  commonly  found  in  the  literature. 

In  Equation  (3.85),  /3kj  have  been  considered  constants.  Scatchard  (1936)  has  shown 
that  /Ikj  is  a  slowly  varying  function  of  the  ionic  strength.  The  interaction  parameters  may  also 
depend  upon  the  temperature  of  the  system.  However,  Edwards  et  al.  found  this  effect  to  be 
small. 

Pitzer  (1973)  has  derived  a  more  accurate  expression  for  the  activity  coefficients.  His 
expression  is  much  too  complicated  to  be  used  for  alkanolamine  solutions.  The  available 
information  on  the  reactions  of  the  alkanolamines  is  insufficient  to  determine  the  specific 
interaction  parameters  in  the  Pitzer’s  expression.  Beutier  and  Renon  (1978)  have  used  Pitzer’s 
expression  to  calculate  the  VLE  in  NH3-CO2-H2O,  NH3-H2S-H2O  and  NH3-SO2-H2O  systems. 
Though  the  reactions  occurring  in  these  systems  are  well-studied,  Beutier  and  Renon  had  to 
make  several  assumptions  about  the  interaction  parameters  to  make  the  expression 
manageable.  The  extended  Debye-Huckel  expression  seems  to  be  the  appropriate  one  for  the 
CO2  /  H2S  /  amine  /  water  system.  Given  the  complexity  of  this  expression,  a  good  computer 
program  for  solving  the  VLE  equations  becomes  a  necessity. 


IV.  METHOD  OF  SOLUTION 


The  problem  of  correlating  the  VLE  data  consists  of  (1)  the  interpretation  of  the 
data  to  obtain  the  required  parameters  and  (2)  the  application  of  the  knowledge  thus  obtained 
to  predict  the  equilibria.  In  both  stages,  one  needs  a  numerical  method  to  relate  the  vapor  and 
the  liquid  compositions. 

The  mathematical  problem  is  to  solve  Equations  (3.11)  to  (3.24).  The  fugacity 
coefficients  in  these  equations  are  obtained  from  the  Peng-Robinson  equation  of  state.  The 
activity  coefficients  are  obtained  from  the  extended  Debye-Huckel  theory  of  electrolyte 
solutions.  Essentially  one  has  to  solve  a  system  of  nonlinear  algebraic  equations. 

The  sequence  of  calculations  is  as  follows:  (1)  The  temperature,  the  normality  of 
amine  and  either  the  acid  gas  partial  pressures  or  the  mole  ratios  in  liquid  are  specified. 

(2)  The  equilibrium  constants,  the  Henry’s  constants,  the  dielectric  constant  of  water  and  the 
binary  interaction  parameters  are  calculated.  (3)  A  preliminary  solution  is  calculated  by  setting 
all  activity  coefficients  and  the  fugacity  coefficients  equal  to  unity.  (4)  The  preliminary  estimate 
is  refined  using  Brown’s  method  for  solving  a  system  of  nonlinear  equations.  (5)  The 
equilibrium  composition  is  determined  from  the  converged  solution. 

A  FORTRAN  program  for  calculating  the  equilibrium  composition  is  presented  in 
Appendix  C.  Brown’s  method  used  in  this  program  is  based  on  a  partial  pivoting  technique 
similar  to  Gaussian  elimination.  The  method  is  at  least  second  order  convergent,  requires  fewer 
function  evaluations  as  compared  to  the  Newton’s  method  and  is  derivative-free.  This  method 
is  especially  attractive  when  the  number  of  variables  is  large  or  when  the  equations  to  be 
solved  are  complicated.  Brown’s  method  is  described  in  Appendix  B. 

To  obtain  a  fast  convergence  with  Brown’s  method,  the  variables  and  the  functions 
have  to  be  “normalized”.  The  variables  are  normalized  if  they  are  all  of  the  same  order  of 
magnitude.  The  molalities  and  the  partial  pressures  in  the  CO2  /  H2S  /  amine  /  water  system 
range  over  several  orders  of  magnitude.  Molalities  as  low  as  10  10  mol/kg-water  and  pressures 
as  high  as  106  Pa  are  common.  To  normalize  these  variables,  the  equations  were  formulated  in 
terms  of  their  logarithms.  For  example,  instead  of  p(CC>2),  In  p(CC>2)  was  treated  as  a  variable. 
Even  though  the  original  variables  range  over  several  orders  of  magnitudes,  their  logarithms 
are  of  the  same  order  of  magnitude.  When  the  logarithms  are  treated  as  variables,  the  partial 
derivatives  in  the  Jacobian  matrix  are  calculated  with  respect  to  the  logarithms.  Instead  of  the 
original  problem,  f(x)  =  0,  we  solve  the  problem  g(y)  =  0;  where  y,  =  In  x,.  The  use  of 
logarithms  has  an  additional  advantage  that  the  solution  remains  feasible  at  each  step. 
Regardless  of  the  magnitude  and  the  sign  of  yi,  Xj  always  remains  positive. 

The  functions  are  normalized  if  the  convergence  criterion  can  be  correctly  expressed 
as 

|f(x)|  <  e 

where  e  is  a  reasonably  small  positive  constant.  For  example,  consider  the  function 
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f(Xi,  X2)  =  Xi  -  x2 
and  the  error  criterion 
|f(xi,  x2)|  <  10~3. 

If  xi  =  100  000  and  x2  =  99  999,  the  error  criterion  is  not  satisfied  even  though  Xi  and  x2 
agree  to  four  digits.  On  the  other  hand,  if  xi  =  10~15  and  x2  =  10~20,  the  error  criterion  is 
satisfied  even  though  xi  is  100  000  times  as  large  as  than  x2.  This  means  that  the  function 
f(xj,  x2)  is  not  normalized.  A  better  way  to  write  this  function  would  be 
g(xi,  x2)  =  1  -  Xi/X2 

In  this  case,  the  error  criterion  would  be  meaningful  for  any  Xi  and  x2. 

Equations  (3.11)  to  (3.24)  may  be  normalized  as  follows: 
fi  =  1  -  7(H+)  7(OH")  m(H+)  m(OH")  /  Kw  xw 
f2  =  1  -  </>( C02)  y(C02)  P/H(C02)  7(C02)  m(C02) 
f3  =  1  -  (  m(H2S)  +  m(HS~)  +  m(S=)  )  /  a(H2S)  mA 
f4  =  1  -  S  (negative  charges)  /  2  (positive  charges) 

The  remaining  equations  may  be  normalized  in  the  same  way. 

After  the  variables  and  the  functions  are  thus  normalized,  the  programming  is 
generally  straightforward,  unless  one  or  more  components  are  absent.  For  example,  if  no 
hydrogen  sulfide  is  present  in  the  system,  m(H2S),  m(HS^  and  m(S=)  are  all  equal  to  zero.  As 
In  0  is  undefined,  the  method  described  above  has  to  be  modified  to  handle  these  cases.  This 
is  done  by  using  two  arrays,  Iv  and  In.  These  arrays  are  defined  as  follows: 

The  array  Iv  contains  the  indices  of  the  variables  to  be  determined.  When  a  variable 
is  known  to  be  zero,  its  index  is  not  included  in  Iv.  If  x2  =  0,  Iv  =  [1,  3,  4,  .  .  .  ].  Similarly, 
the  array  Ie  contains  the  indices  of  the  equations  to  be  solved.  These  arrays  can  be  calculated 
if  the  species  present  in  the  system  are  known.  If  all  components  are  present  in  the  system, 

Iv(i)  =  Ie(  i)  =  i. 

If  some  components  are  absent  from  the  system,  a  smaller  set  of  equations  is  solved. 
The  equations  solved  are: 

fk(x)  =  0  k  =  IE(1),  Ie(2),  .  .  . 

These  equations  are  solved  in  terms  of  the  variables  x,  where 

j  =  Iv(l),  Iv<2),  .  .  . 

Using  this  technique,  we  avoid  taking  logarithms  of  quantities  which  are  zero,  while 
normalizing  the  non-zero  variables. 

The  indexing  also  permits  us  to  use  the  same  program  to  calculate  the  partial 
pressures  when  the  loadings  are  known  and  to  calculate  the  loadings  when  the  partial 
pressures  are  known.  If  the  partial  pressure  of  C02  is  known,  the  index  of  the  variable 
denoting  p(C02)  is  not  included  in  Iv.  At  the  same  time,  the  index  of  the  variable  denoting 
a(C02)  is  included  in  Iv.  This  has  the  effect  of  treating  a(C02)  as  an  unknown  and  p(C02)  as 
a  known  variable.  When  a(C02)  is  known  and  p(C02)  is  unknown,  the  index  of  a(C02)  is 
excluded  from  Iv  and  that  of  p(C02)  is  included  in  Iv. 

In  summary,  the  computer  program  has  the  following  features: 

1.  The  algorithm  is  fast  and  efficient. 

2.  The  program  is  able  to  calculate  the  partial  pressures  if  the  mole  ratios  in  the  liquid 
phase  are  known  and  vice  versa. 
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3.  Although  the  variables  are  normalized  by  taking  logarithms,  the  program  can  handle 
variables  which  are  zero.  Redundant  equations  are  also  avoided. 

If  another  method  is  to  be  used  to  calculate  the  fugacity  coefficients  and  the  activity 
coefficients,  the  program  can  be  easily  modified.  Only  the  routines  that  calculate  these 
coefficients  have  to  be  rewritten. 

With  this  program,  we  can  predict  the  vapor-liquid  equuilibria  if  the  equilibrium 
constants,  the  Henry’s  constants  and  the  interaction  parameters  are  known.  There  are  no  a 
priori  methods  to  determine  the  specific  interaction  parameters.  These  parameters  have  to  be 
deduced  from  the  experimental  data. 

The  specific  interaction  parameters  show  a  weak  dependence  on  the  ionic  strength 
(Scatchard,  1936)  and  the  temperature  (Edwards  et  a!.,  1975).  We  assume  these  parameters  to 
be  independent  of  the  temperature  and  the  ionic  strength. 

The  specific  interaction  parameters  represent  the  short  range  interactions  (van  der 
Waals  forces).  In  moderately  charged  solutions,  the  similarly  charged  ions  do  not  come  close 
enough  for  the  van  der  Waals  forces  to  become  important.  This  is  called  the  Br^nsted 
principle  (Br^»nsted,  1922  and  1923).  If  i  and  j  are  similarly  charged  (both  positive  or  both 
negative)  ions,  the  Br^nsted  principle  implies  that  (3, y  =  0. 

The  interaction  parameters  may  be  divided  into  three  groups:  (1)  the 
molecule-molecule  interaction  parameters,  (2)  the  molecule-ion  interaction  parameters  and 
(3)  the  ion-ion  interaction  parameters. 

Edwards  et  al.  have  determined  /Ta,  the  molecule-molecule  interaction  parameters  by 
fitting  Equation  (2.85)  to  single-solute  systems.  As  p  —  0, 

In  7a  2  /Saa  UTa 

Thus  (2  /Ta)  is  the  limiting  slope  of  the  plot  of  In  y a  against  m*.  However,  as  the 
concentration  of  the  solute  tends  to  zero,  the  degree  of  dissociation  tends  to  infinity.  At  zero 
ionic  strength,  the  electrolyte  will  be  fully  dissociated.  Edwards  et  al.  assumed  that  the 
molecule-ion  interaction  parameters  in  the  single  solute  systems  are  equal  to  zero.  Without  this 
assumption,  /3aa  can  not  be  determined  by  Edwards’  procedure.  The  molecule-molecule 
interaction  parameters  for  CO2,  H2S  and  RR'NH  can  be  determined  from  the  data  of 
Houghton  et  al.  (1957),  Clarke  and  Glew  (1970)  and  Sivertz  et  al.  (1940).  The 
molecule-molecule  interaction  parameters  between  the  unlike  species  may  be  estimated  from  the 
combination  rule  proposed  by  Edwards  et  al.: 

/Sab  =  0.5  (  /?a  +  (3b  ) 

As  no  data  are  available  for  DIPA,  DGA  and  MDEA,  it  is  proposed  that 
/Saa(DGA)  =  /Saa(MEA) 

/Saa(DIPA)  =  /Saa(DEA) 

/Saa(MDEA)  =  /Saa(TEA) 

Bromley  (1972)  has  formulated  a  procedure  for  estimating  the  binary  ion-ion 
interaction  parameters.  He  postulated  that  the  interaction  parameter  /S,j  may  be  expressed  as  a 
sum  of  individual  ion  contributions,  (3t  and  f3j. 

Pa  =  Pi  +  P, 

Further,  Bromley  observed  a  correlation  between  (3i  and  the  standard  partial  molar  entropy  of 
the  ion,  z,Si°,  at  25°C.  If  the  standard  entropy  of  the  ions  are  known,  /3,  and  pti  may  be 
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estimated  from  Bromley’s  plot. 

The  procedure  works  well  for  CO2  /  H2O  and  H2S  /  H2O  systems.  The  entropies  of 
ethanolammonium  and  the  carbamate  ions  have  not  been  measured.  Bromley’s  procedure  can 
not  be  applied  for  these  ions.  Bromley  (1973)  has  also  proposed  a  more  accurate  method  for 
estimating  the  ion-ion  interaction  parameters.  In  addition  to  standard  entropies,  this  method 
requires  the  knowledge  of  the  proton  affinities  of  the  ions. 

Beutier  and  Renon  have  found  out  that  a  unique  correlation  between  /?,  and  z,Si° 
does  not  exist  for  all  ions  in  solution.  They  observed  different  correlations  for  cations, 
halogenide  ions,  oxygenated  polyatomic  ions  and  proton  acceptors.  Even  if  the  entropies  of  the 
ethanolammonium  and  the  carbamate  ions  could  be  estimated,  Bromley’s  plot  could  not  be 
used  with  confidence. 

The  molecule-ion  (or  salting-out)  parameters  may  also  be  correlated  in  this  way. 
Edwards  et  al.  assumed  that  pm\  depends  linearly  on  p,  or  ZjS,°.  Long  and  McDevit  (1952) 
have  measured  the  salting-out  parameters  for  carbon  dioxide  and  hydrogen  sulfide.  Their  data 
can  be  used  to  construct  a  plot  of  pmi  against  ZiSi0.  This  method  can  be  used  for  carbon 
dioxide  /  hydrogen  sulfide  /  water  system.  The  salting-out  parameters  for  amine-containing 
systems  have  not  been  measured.  These  parameters  have  to  be  obtained  from  least-squares 
fitting  of  the  experimental  data. 

Let  B  =  [jSi,  Pi,  •  •  •  ,  /?n]  be  the  array  of  the  interaction  parameters  and  E  =  [ei, 

62,  •  •  •  ,  Cm]  be  the  array  of  the  errors  in  predictions.  We  wish  to  find  the  interaction 
parameters  B*  that  minimize  s  =  |E|2.  When  s  is  minimum, 

3s  /  dp,  =  0,  i  =  1,  2,  .  .  .  ,N 

These  N  equations  can,  in  principle,  be  solved  for  the  N  unknown  P's.  However,  such  a 
procedure  is  not  feasible  because  of  the  excessive  computational  effort  involved.  To  simplify 
the  calculations,  the  following  assumptions  may  be  made. 

1.  The  interaction  parameters  for  the  species  which  are  present  in  very  small  concentrations 
can  be  neglected. 

1.  The  effect  of  the  P' s  on  the  errors  in  prediction  is  assumed  to  be  additive. 

When  a  certain  chemical  species  is  present  in  a  very  small  concentration,  it  has  little  effect  on 
the  equilibria.  The  concentration  of  that  species  may  be  neglected  from  the  balance  equations. 
A  reduced  set  of  equations  can  be  solved  in  terms  of  a  smaller  set  of  variables.  For  example, 
if  S=  is  present  in  a  very  small  amount,  its  concentration  can  be  ignored  from  the  charge 
balance  and  the  sulfur  balance  equations.  The  equations  (except  the  equilibrium  relation  for 
the  dissociation  of  the  bisulfide  ion)  can  be  solved  for  the  variables  (excluding  S').  Secondly, 
the  term  /^nij  in  the  extended  Debye-Huckel  theory  may  be  neglected  if  m,  is  small. 

The  additivity  assumption  implies  that  the  total  deviation  in  predictions  may  be 
expressed  as  a  sum  of  the  individual  deviations  due  to  each  of  the  interaction  parameters.  Let 
P°  be  the  predictions  using  B  =  [0,  0,  .  .  .  ,  0]  and  let  P  be  the  measured  values.  Therefore, 

E°  =  p  -  p°  is  the  total  error.  The  effect  of  the  interaction  parameters  can  be  found  by  setting 
j-th  parameter  to  a  small  number,  say  Ap,.  Let  the  new  predictions  be  P(A/3j).  This  means  that 
setting  the  j-th  parameter  to  A/Sj  has  caused  a  deviation  of  [P(  A/3j)-P°]  in  the  predictions.  If 
P's  are  small,  this  deviation  is  proportional  to  the  magnitude  of  P ;  i.e. 

[P(ft)  -  p °]  /  ft  =  [P(AA)  -  P°]  /  Aft 


' 
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P08j)  -  P°  =  08j  /  A/3j)  [P(A/3j)  -  P0] 

Define  C,  =  [P(A/3;)  -  P°]/A^j.  Since  the  total  deviation  is  the  sum  of  the  individual 
deviations, 

CT  =  2  ft  Q 
This  can  be  written  as 

CT  =  X  B 

where  X  =  [Ci,  C2,  .  .  .  ,  Cn].  We  wish  to  find  B  so  that  Ct  is  the  best  estimate  to  E°,  i.e. 
we  wish  to  minimize  the  squared  error  q  =  |E°  -  Ct|2. 

q  =  (E°  -  Ct)'-(E°  -  CT) 

Substituting  Ct  =  X  B  and  expanding, 

q  =  E°'E°  -  B'X'E0  -  E°'XB  +  B'X'XB 
dq/dfii  =  -2  l'X'E0  +  2  l'X'XB 
When  q  is  minimum,  9q/9ft  =  0. 

..  X  XB  =  X  E° 

B  =  (X'Xr’X'E0 

All  arrays  required  in  this  method  can  be  calculated  by  using  the  the  computer 
program  described  above.  The  method  can  be  extended  to  determine  the  equilibrium  constants 
outside  their  range  of  measurement.  Define  K  =  Ko  exp  e.  The  exponential  function  is  used  to 
ensure  that  K  is  always  positive.  Given  an  initial  estimate  of  K,  e  can  be  determined  by  the 
least  squares  fitting  as  described  above. 

This  method  was  used  to  correlate  the  equilibrium  solubility  of  carbon  dioxide, 
hydrogen  sulfide  and  their  mixtures  in  MEA,  DEA  and  DIPA  solutions.  The  results  are 
presented  in  the  next  chapter. 


V.  RESULTS  AND  DISCUSSION 


The  procedure  described  in  Chapter  IV  was  applied  to  the  CO2  /  H2S  /  amine  / 
water  system.  The  temperature,  amine  normality  and  the  acid  gas  loadings  were  provided  as 
the  input  to  the  computer  program  and  the  predicted  partial  pressures  were  compared  to  the 
experimental  data.  The  error  in  predictions  was  taken  to  be  In  P0bS  -  In  Ppred.  The  interaction 
parameters  were  determined  by  minimizing  the  sum  of  the  squared  errors. 

In  Brown’s  method,  the  iterations  were  stopped  when  two  successive  estimates  of  the 
solution  agreed  to  within  eight  significant  digits  or  when  all  functions  were  less  than  10~15  in 
absolute  value.  These  criteria  appear  to  be  unnecessarily  stringent  at  first  sight.  Brown’s 
method  possesses  second  order  convergence  properties,  i.e.  the  number  of  significant  digits 
approximately  doubles  in  each  step  if  the  estimate  is  close  enough  to  the  solution.  If  it  takes 
five  iterations  to  obtain  a  solution  correct  to  two  significant  places,  it  will  take  just  two  more 
iterations  to  obtain  an  estimate  correct  to  eight  significant  places.  Thus  the  extra  effort 
required  to  increase  the  number  of  significant  digits  is  quite  small.  Secondly,  making  the  error 
criterion  so  stringent  ensures  that  no  errors  are  introduced  by  the  numerical  method  itself. 

In  the  CO2  /  H2S  /  amine  /  water  system,  the  only  species  which  are  present  in 
significant  concentration  are  the  free  amine  and  the  ethanolammonium,  carbamate,  bicarbonate 
and  bisulfide  ions.  This  means  that  the  interaction  parameters  for  the  other  species  (C02,  ITS, 
OH“,  H+,  S=  and  CO3)  will  have  little  effect  on  the  observed  equilibria.  At  very  high  loadings, 
the  concentrations  of  the  free  CO2  and  H2S  are  also  high. 

The  important  interaction  parameters  are  /3(RR'NH-RR'NH), 

^(RR'NH-RR'NCOOl,  /3(RR'NH-HC03T  /3(RR'NH-HS^,  £(RR'NH2+-RR'NCOCr), 
£(RR'NH2+-HC03~)  and  /3(RR'NH2+-HS~).  The  first  of  these  can  be  determined  from  the 
correlation  of  Atwood  et  al..  It  should  be  noted  that  there  are  no  mixture  parameters  in  this 
list.  The  parameters  for  the  C02  /  amine  /  water  system  and  the  H2S  /  amine  /  water  system 
are  sufficient  to  predict  the  the  equilibria  in  the  CO2  /  ITS  /  amine  /  water  system.  The 
required  parameters  may  be  obtained  either  from  the  single  solute  data  or  the  mixture  data. 
More  computations  are  required  if  the  mixture  data  are  used  to  calculate  the  parameters. 


1.  Solubility  of  acid  gases  in  MEA  solutions 

The  partial  pressures  of  ITS  over  2.5  N  and  5.0  N  MEA  solutions  at  25,  40,  60,  80, 
100  and  120°C  have  been  reported  by  Lee  et  al.,  1976a.  Their  data  were  used  to  determine 
/}(RR'NH-HST  /3(RR'NIT+-HS^)  and  the  equilibrium  constant  for  the  dissociation  of  MEA  at 
100°C.  Figures  2  and  3  show  comparisons  between  the  predictions  and  the  experimental  data. 
The  agreement  is  good  for  both  2.5  N  and  5.0  N  MEA  solutions  and  at  all  temperatures. 

Lee  et  al.,  1976b  have  measured  the  solubility  of  C02  in  1.0,  2.5,  3.75  and  5.0  N 
MEA  solutions  at  25,  40,  60,  80,  100  and  120°C.  The  equilibrium  constant  for  the  carbamate 
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H2S  PART  I RL  PRESSURE  ,  KPfl 
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MOLE  RRTIO  IN  LIQUID,  H2S/MER 


Figure  2. 


SOLUBILITY  OF  H2S  IN  2-5  N  MEfl  SOLUTIONS 


H2S  PARTIAL  PRESSURE  ,  KPA 
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Figure  3. 


SOLUBILITY  OF  H2S  IN  5-0  N  MEfl  SOLUTIONS 


C02  PRRTIRL  PRESSURE,  KPR 
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Figure  4. 


SOLUBILITY  OF  C02  IN  1.0  N  HER  SOLUTIONS 


C02  PARTIAL  PRESSURE  ,  KP A 
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MOLE  RATIO  IN  LIQUID,  C02/MEA 


Figure  5. 


SOLUBILITY  OF  C02  IN  2.5  N  MER  SOLUTIONS 


C02  PARTIAL  PRESSURE  ,  KP A 
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MOLE  RATIO  IN  LIQUID,  C02/MEA 


Figure  6. 


SOLUBILITY  OF  C02  IN  3-75  N  MEfl  SOLUTIONS 


C02  PARTIAL  PRESSURE  ,  KP A 
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MOLE  RATIO  IN  LIQUID,  C02/MEA 


Figure  7. 

SOLUBILITY  OF  C02  IN  5-0  N  MEfl  SOLUTIONS 


H2S  PARTIAL  PRESSURE,  KPA 
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MOLE  RATIO  IN  LIQUID,  H2S/MEA 


Figure  8. 


EFFECT  OF  C02 
IN  2.5  N  MER 


ON  THE  SOLUBILITY  OF  H2S 
SOLUTIONS  RT  40  CELSIUS 


C02  PARTIAL  PRESSURE,  KP A 
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Figure  9. 


OF  H2S  ON  THE  SOLUBILITY  OF  C02 
N  ME A  SOLUTIONS  AT  40  CELSIUS 


H2S  PART  I RL  PRESSURE  ,  KPR 
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MOLE  RRTIO  IN  LIQUID,  H2S/ME- 


Figure  10. 

EFFECT  OF  C02  ON  THE  SOLUBILITY  OF  -23 
IN  5.0  N  MER  SOLUTIONS  RT  40  CELSIUS 


C02  PARTIAL  PRESSURE  »  KP A 
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EFFECT  OF  H2S  ON  THE  SOLUBILITY  OF  C02 
IN  5.0  N  MER  SOLUTIONS  RT  40  CELSIUS 


H2S  PARTIAL  PRESSURE  ,  KP A 
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MOLE  RATIO  IN  LIQUID  ,  H2S/MEA 


Figure  12. 

EFFECT  OF  C02  ON  THE  SOLUBILITY  OF  H2S 
IN  2.5  N  MER  SOLUTIONS  RT  100  CELSIUS 


C02  PARTIAL  PRESSURE  ,  KPA 
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MOLE  RATIO  IN  LIQUID  ,  C02/MEA 


Figure  13. 


EFFECT  OF  H2S  ON  THE  SOLUBILITY  OF  C02 
IN  2-5  N  MER  SOLUTIONS  AT  100  CELSIUS 


H2S  PARTIAL  PRESSURE  ,  KPA 
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MOLE  RATIO  IN  LIQUID  ,  H2S/MEA 

Figure  14. 

EFFECT  OF  C02  ON  THE  SOLUBILITY  OF  H2S 
IN  5.0  N  MER  SOLUTIONS  RT  100  CELSIUS 


C02  PRRTIRL  PRESSURE,  KPR 
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MOLE  RRTIO  IN  LIQUID,  C02/MER 


Figure  15. 


EFFECT  OF  H2S  ON  THE  SOLUBILITY  OF  C02 
IN  5.0  N  MEfl  SOLUTIONS  HT  100  CELSIUS 
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formation  at  100°C,  as  well  as  the  interaction  parameters  for  the  CO2  /  amine  /  water  system 
were  determined  from  their  data.  Figures  4,  5,  6  and  7  compare  the  predictions  to  the 
experimental  data.  The  agreement  is  good  at  all  normalities  and  temperatures. 

The  constants  obtained  from  the  single  solute  data  were  used  to  predict  the  equilibria 
for  mixtures  of  carbon  dioxide  and  hydrogen  sulfide.  Comparisons  with  the  experimental  data 
of  Lee  et  al.,  1976c  (2.5  N  MEA  solutions)  and  Lee  et  al.,  1975  (5.0  N  MEA  solutions)  are 
shown  in  Figures  8  to  14.  The  agreement  with  experiment  is  much  better  at  100°C  than  at 
40°C,  probably  because  the  carbamate  formation  is  unimportant  at  the  higher  temperature.  At 
40°C,  the  agreement  gets  worse  as  the  loading  of  the  other  acid  gas  increases.  This  is  only  to 
be  expected  because  no  mixture  data  have  been  used  to  determine  the  constants.  It  also 
suggests  that  the  interaction  parameters  such  as  /3(C02-HS~)  and  /3(H2S-HC03^)  become 
important  at  high  loadings.  On  the  whole,  the  agreement  is  much  better  than  that  of  the 
previous  models.  Generally,  the  predictions  tend  to  bracket  the  experimental  data.  At  loadings 
exceeding  unity,  however,  the  predicted  partial  pressures  are  usually  lower  than  the  the 
experimental  values. 

The  most  likely  direction  for  improving  the  predictions  seems  to  be  the  inclusion  of 
parameters  involving  CO2  and  H2S.  This  would  require  the  equilibrium  data  for  the  CO2  / 

H2S  /  water  system.  Unfortunately,  there  are  very  few  data  available  for  this  system.  It  would 
be  very  helpful  to  have  reliable  experimental  data  for  the  carbamate  formation. 


2.  Solubility  of  acid  gases  in  DIP  A  solutions 

Isaacs  et  al.  (1977a  and  1977b)  have  reported  the  solubility  of  carbon  dioxide, 
hydrogen  sulfide  and  their  mixtures  in  2.5  N  DIPA  solutions  at  40  and  100°C.  Their  data  for 
the  pure  components  as  well  as  mixtures  were  used  to  deduce  the  binary  interaction 
parameters  and  the  equilibrium  constants.  These  are  presented  in  Table  2. 

The  fitted  equilibrium  constants  for  DIPA  are  of  the  same  order  magnitude  as  the 
measured  equilibrium  constants  for  MEA  and  DEA.  Figures  16,  17,  18  and  19  compare  the 
predictions  made  by  using  these  constants  to  experimental  data.  As  compared  to  MEA,  there 
is  a  noticeable  improvement  in  accuracy,  particularly  for  mixtures.  The  predictions  at  very  low 
mole  ratios  in  the  liquid  phase  (less  than  0.1)  are  not  as  good  as  those  at  moderate  or  high 
loadings.  This  fact  can  not  be  considered  significant  because  much  of  the  data  for  low 
loadings  were  obtained  by  extrapolation  and  not  by  direct  measurement. 


3.  Solubility  of  acid  gases  in  DEA  solutions 

Lee  et  al.  have  measured  the  solubility  of  mixtures  of  carbon  dioxide  and  hydrogen 
sulfide  in  DEA  solutions.  Their  data  for  the  solubility  of  carbon  dioxide  and  hydrogen  sulfide 
were  used  to  obtain  the  binary  interaction  parameters  and  the  equilibrium  constants.  These 
constants  are  presented  in  Table  2.  Figures  20  to  27  show  the  predictions  and  the  experimental 
solubility  data  for  2.0  and  3.5  N  DEA  solutions  at  50  and  100°C. 

The  predictions  of  the  solubility  of  the  carbon  dioxide  in  3.5  N  DEA  solutions  are 


H2S  PRRTIRL  PRESSURE  ,  KPR 
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MOLE  RATIO  IN  LIQUID  ,  H2S/DIPR 


Figure  16. 


EFFECT  OF  C02  ON  THE  SOLUBILITY  OF  H2S 
IN  2.5  N  DIPfl  SOLUTIONS  AT  40  CELSIUS 


C02  PARTIAL  PRESSURE  ,  KPA 
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MOLE  RATIO  IN  LIQUID  ,  C02/DIPA 

Figure  17. 

EFFECT  OF  H2S  ON  THE  SOLUBILITY  OF  C02 
IN  2.5  N  D I  PR  SOLUTIONS  AT  40  CELSIUS 


H2S  PARTIAL  PRESSURE,  KPA 
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MOLE  RATIO  IN  LIQUID,  H2S/DIPA 


Figure  18. 


EFFECT  OF  C02  ON  THE  SOLUBILITY  OF  H2S 
IN  2.5  N  DIPfl  SOLUTIONS  AT  100  CELSIUS 


C02  PARTIAL  PRESSURE,  KPR 
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MOLE  RATIO  IN  LIQUID,  C02/DIPA 


Figure  19. 


EFFECT  OF  H2S  ON  THE  SOLUBILITY  OF  C02 
IN  2.5  N  D I  PR  SOLUTIONS  RT  100  CELSIUS 


H2S  PARTIAL  PRESSURE  ,  KPA 
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MOLE  RATIO 


IN  LIQUID,  H2S/DEA 


Figure  20. 


EFFECT  OF  C02  ON  THE  SOLUBILITY  OF  H2S 
IN  2.0  N  DEfl  SOLUTIONS  AT  50  CELSIUS 


C02  PRRTIRL  PRESSURE  ,  KPR 
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MOLE  RRTIO  IN  LIQUID  ,  C02/DER 


Figure  21. 


EFFECT  OF  H2S  ON  THE  SOLUBILITY  OF  002 
IN  2.0  N  DEfl  SOLUTIONS  AT  50  CELSIUS 


H2S  PARTIAL  PRESSURE  ,  KPA 
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EFFECT  OF  C02  ON  THE  SOLUBILITY  OF  H2S 
IN  3.5  N  DEfl  SOLUTIONS  F)T  50  CELSIUS 


C02  PRRTIflL  PRESSURE,  KPR 
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MOLE  RATIO  IN  LIQUID,  002/DEfl 


Figure  23. 


EFFECT  OF  H2S  ON  THE  SOLUBILITY  OF  002 
IN  3-5  N  DEfl  SOLUTIONS  AT  50  CELSIUS 


H2S  PARTIAL  PRESSURE  ,  KPA 
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MOLE  RATIO  IN  LIQUID  ,  H2S/DEA 


Figure  24. 


EFFECT  OF  C02  ON  THE  SOLUBILITY  OF  H2S 
IN  2.0  N  QEfl  SOLUTIONS  AT  100  CELSIUS 


C02  PARTIAL  PRESSURE  ,  KPA 
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MOLE  RATIO  IN  LIQUID,  C02/DEA 


Figure  25. 


EFFECT  OF  H2S  ON  THE  SOLUBILITY  OF  C02 
IN  2.0  N  QEfi  SOLUTIONS  HT  100  CELSIUS 


H2S  PARTIAL  PRESSURE  ,  KPA 
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MOLE  RATIO  IN  LIQUID,  H2S/DEA 

Figure  26. 

EFFECT  OF  C02  ON  THE  SOLUBILITY  OF  H2S 
IN  3.5  N  DEfl  SOLUTIONS  AT  100  CELSIUS 


C02  PARTIRL  PRESSURE  ,  KPA 
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MOLE  RATIO  IN  LIQUID,  C02/DEA 


Figure  27. 


EFFECT  OF  H2S  ON  THE  SOLUBILITY  OF  C02 
IN  3.5  N  DEfl  SOLUTIONS  AT  100  CELSIUS 
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poor.  Other  predictions  are  in  close  agreement  with  the  experimental  data.  As  with  MEA,  the 
agreement  gets  poorer  as  the  mole  ratio  of  the  other  acid  gas  increases.  It  is  likely  that  the 
agreement  will  be  improved  if  the  mixture  data  is  used  to  determine  the  interaction 
parameters. 
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Table  2.  Specific  interaction  parameters 


Monoethanolamine 


0(RR'NH-HS>-O.O375, 
/3(RR'NH2+-HS>0.005, 
0(RR'NH-HCO3>-O.O245, 
/3(RR'NH2+-HCO3~)=-0.0375, 
/?(RR  'NH-RR'NCOO~)= -0.063, 
/3(RR'NH2+-RR'NCOO~)=0.006, 


Diethanolamine 


/3(RR'NH-HS”)=0. 029845, 
/3(RR'NH2+-HS~)=0. 00482, 
0(RR'NH-HCO3~)=O.OO534, 
/3(RR  'NH2+-HC03>=-0.04 1 , 
0(RR'NH-RR'NCOO>=O.26O1, 
0(RR'NH2+-RR'NCOO1=-O.O4288, 


Diisopropanolamine 


0(RR'NH-HS~)=O.O921, 
£(RR'NH2+-HS>0.0905, 
0(RR'NH-HCO3~)=O.3O9, 
/3(RR'NH2+-HC03>0.103, 
/3(RR'NH-RR'NCOO”)= -0.00294, 
yS(  RR  'NH2+-RR  'NC001=0. 1098, 


' 


VI.  CONCLUSIONS 


The  solubility  model  presented  here  is  consistent  with  the  thermodynamic  criteria  for 
the  equilibrium  in  multicomponent,  multiphase  systems.  The  binary  interaction  parameters 
required  in  the  model  are  easily  determined  from  the  experimental  data.  The  equilibrium 
constants  can  also  be  determined  from  the  solubility  data.  The  equilibrium  data  for  the  pure 
components  or  for  the  mixtures  may  be  used  to  determine  the  interaction  parameters.  The 
computer  program  developed  in  this  work  is  efficient  and  general.  It  can  be  used  to  compute 
the  equilibrium  composition  in  any  multicomponent  system  if  the  equilibrium  constants, 

Henry’s  constants  and  the  interaction  parameters  are  known. 

The  predictions  made  using  this  model  are  better  than  those  of  the  Kent  and 
Eisenberg  and  the  Klyamer  et  al.  models.  The  agreement  for  the  mixtures  is  improved  if  the 
mixture  data  are  used  to  calculate  the  interaction  parameters. 

The  model  uses  the  dielectric  constant  of  pure  water  in  the  Debye-Huckel  expression 
where  the  dielectric  constant  of  the  solution  is  required.  For  this  assumption  to  be  valid,  the 
concentrations  of  the  solute  species  must  remain  low.  The  activity  coefficient  of  water  has  been 
assumed  to  be  unity.  This  requires  the  mole  fraction  of  water  in  the  liquid  phase  to  be  close 
to  unity.  The  validity  of  the  model  also  depends  upon  the  assumed  reaction  mechanism.  For 
example,  if  the  pH  is  more  than  11,  the  hydroxyl  group  of  the  alkanolamines  also  reacts  with 
the  acid  gases. 

Many  of  the  equilibrium  constants  required  in  the  model  have  been  measured  only  in 
a  narrow  range  of  temperatures.  Experimental  measurements  of  these  constants  upto  120°C  are 
certainly  desirable.  Solubility  data  for  the  CO2  /  H2S  /  water  system  and  for  the  C02  /  H2S 
/  amine  /  water  system  are  also  needed  to  verify  the  validity  of  the  model  for  these  systems. 
More  binary  interaction  parameters  should  be  included  in  the  model  to  improve  the 
predictions  for  the  high  loadings. 
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APPENDIX  A 


Klyamer  et  al.  Model 


The  Klyamer  et  al.  model  is  based  on  the  assumption  that  the  activity  coefficients  of 
all  ions  in  the  liquid  phase  are  equal.  The  expressions  for  the  partial  pressures  are  derived 
below. 

Reactions: 

RR'NH  +  H20  =  RR'NH2+  +  OH~ 

H2S  =  H+  +  HS~ 

H20  =  H+  +  OH' 

HS"  =  H+  +  S= 

2  RR'NH  +  C02  =  RR'NH2+  +  RR'NCOO' 

H20  +  C02  =  H+  +  HC03' 

HCO3'  =  H+  +  C03= 

Conservation  equations: 

[RR'NH21  +  [Hi  -  [HCO3I  +  [RR'NCOOl  +  [OH1  +  [HS1 


+  2  [C03=]  +  2  [S=]  —  (A. 62) 

mA  =  [RR'NH]  +  [RR'NH21  +  [RR'NCOOl  —  (A.63) 

mA  x(C02)  -  [CO2]  +  [HCO3I  +  [RR'NCOOl  +  [C03=]  —  (A.64) 

mA  x(H2S)  =  [H2S]  +  [HS1  +  [S=]  —  (A.65) 

K,  =  y2  [RR'NH21  [OH1  /  a  a2  [RR'NH]  —  (A.66) 

Kw  -  y2  [Hi  [OH1  /  a  -  (A.67) 

Kic  -  y2  [Hi  [HS1  /  [H2S]  -  (A.68) 

K2C  =  y  [Hi  [S=]  /  [HS1  -  (A. 69) 

H'(H2S)  =  [H2S]  /  p(H2S)  —  (A. 70) 

Km  =  y2  [RR'NH21  [RR'NCOOl  /  a2  a2  [RR'NH]2  p(C02)  —  (A.71) 

K,y  =  y2  [Hi  [HCO3I  /  a  [CO2]  -  (A. 72) 

K2Y  =  y  [Hi  [OV]  /  [HCO3I  —  (A. 73) 

H'(C02)  =  [CO2]  /  p(C02)  —  (A. 74) 


Derivation  of  expressions  for  calculating  p(COi)  and  pCHiS) 

Define 

A  =  m  x(H2S)  -  H'(H2S)  p(H2S) 

B  =  m  x(C02)  -  H'(C02)  p(C02) 

and 

z  =  [RR'NCOOl 

Neglecting  [Hi,  [OH1,  [CO3I  and  [S=]  in  the  conservation  equations,  and  using  the  above 
notation, 

[HS1  =  m  x(H2S)  -  [H2S]  =  m  x(H2S)  -  H'(H2S)  p(H2S)  =  A 

[HCO3I  =  m  x(C02)  -  [C02]  -  [RR'NCOOl  =  m  x(C02)  -  H'(C02)  p(C02)  -  z  =  B  -  z 
[RRNH2I  =  [HCO3I  +  [RR'NCOOl  +  [HS1  =  B-  z+  z  +  A=  A+  B 


. 
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[RR'NH]  =  m  -  [RR'NHzl  -  [RR'NCOOl  =  m  -  z  -  A  -  B 
Multiply  Equations  (A. 5),  (A. 7),  (A. 9)  and  divide  by  Equation  (A. 6).  On  rearrangement,  we 
get 

p(H2S)  =  (Kw  /  K,KicH'(H2S))  -(y2  /  aa)-A-(A  +  B)  /  (m  -  z  -  A  -  B)  —  (A. 75) 

Rearranging  Equation  (A.  10), 

p(C02)  =  (y2  /  a2a2KM)-z  -(A  +  B)  /  (m  -  z  -  A  -  B)2  —  (A. 76) 

Multiply  Equations  (A. 5),  (A.  11),  (A.  13)  and  divide  by  Equations  (A. 6),  (A.  10)  to  get 

K,K,yH,(C02)  /  KwKm  =  (B  -  z)<m  -  z  -  A  -  B)  /  az  —  (A. 77) 

After  some  manipulations, 

A  =  m  -  z  -  B  -  (KiKiYH'(C02)  /  KWKM)  -z  /  a(B  -  z) 

When  [C02]  <1  [HC031  +  [RR'NCOOl  and  [H2S]  <1  [HS1, 

A  =  m  x(H2S) 

B  =  m  x(C02) 

z  =  B  -  [HCOal  =  m  x(C02)  -  [HC031 
Substitute  in  Equations  (A.  14)  and  (A.  15)  to  obtain 

p(H2S)  =  (Kw  /  K,K1CH'(H2S))  iy2  /  aa)  Xmx(H2S)(x(H2S)  +  x(C02)) 

/  (1  -  x(H2S)  -  2x(C02)  +  [HC031  /  m)  —  (A. 78) 

and 

p(C02)  =  (72  /  KMaV)  X(x(H2S)  +  x(C02))(x(C02)  -  [HCO3I  /  m) 

/  (1  -  x(C02)  -  2x(H2S)  +  [HCO3I  /  m)2  —  (A. 79) 

Substitution  in  Equation  (A.  16)  gives 

mx(H2S)  =  m  -  (mx(C02)  -  [HCO3D  -  mx(C02) 

-  (K,  /  a)(mx(C02)  -  [HCO3D  /  [HCO3I 
where 

K,  =  K,K,yH'(C02)  /  KwKm 
Multiply  the  equation  by  [HCO3I  /  m2  and  rearrange. 

([HCOjI  /  m)2  +  ([HCO3I  /  m)(l  -  x(H2S)  -  2x(C02)  +  K,  /  am) 

-  Kix(C02)  /  am  =  0 

•••  [HCO3I  /  m  =  -  0.5(1  -  x(H2S)  -  2x(C02)  +  K,  /  am) 

±  (  [0.5(1  -  x(H2S)  -  2x(C02)  +  Ki  /  am)]2  +  K,x(C02)  /  am  f 
To  determine  the  sign  of  the  discriminant,  consider  the  case  where  x(C02)  =  0,  x(H2S)  =  0 
and  hence  [HCO3I  =  0-  Negative  sign  of  the  discriminant  leads  to  an  absurd  result,  viz.  1 
+  Ki  /  am  =  0.  As  Ki  /  am  is  positive, 

[HCO3I  /  m  =  (  [0.5(1  -  x(H2S)  -  2x(C02)  +  K,  /  am)]2  +  K,x(C02)  /  am  f5  -  0.5(1 

-  x(H2S)  -  2x(C02)  +  K,  /  am)  —  (A. 80) 

When  [HCO3I  [RR'NCOOl,  Equations  (A.  17)  and  (A.  18)  reduce  to 

p(H2S)  =  (Kw  /  KiKlcH'(H2S))  -(72) 

/  aa)-mx(H2S)(x(H2S)  +  x(C02))  /  (1  -  x(H2S)  -  2x(C02))  —  (A.81) 

p(C02)  =  (72  /  KMa2a2) 

Xx(C02)(x(H2S)  +  x(C02))  /  (1  -  x(H2S)  -  2x(C02))2  —  (A.82) 

When  x(C02)  =  0,  z  =  0,  B  =  0  and  Equation  (A.  14)  becomes 
p(H2S)  =  (Kw  /  K,K1cH'(H2S))  <72  /  aa)-A2  /  (m  -  A) 

Define  K"  =  (Kw  /  K.K.c)-^2  /  aa) 
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•••  p(H2S)  =  (K"  /  H'(H2S))A2  /  (m  -  A) 

Substitute  A  =  mx(H2S)  -  H'(H2S)p(H2S),  expand  and  rearrange  to  get 
w2  +  w(l  -  x(H2S)  +  2K"x(H2S))  =  K"(l  -  K")x2(H2S) 
where  w  =  (1  -  K")p(H2S)H'(H2S)  /  m.  The  solution  to  this  equation  is 

w  +  0.5(1  -  x(H2S)  +  2K"x(H2S))  =  ±  (  [0.5(1  -  x(H2S))]2  +  K"x(H2S)  )05 
The  sign  of  the  discriminant  has  to  be  positive.  (Consider  the  case  x(H2S)  =  0,  p(H2S)  =  0.) 

•••  w  =  (  [0.5(1  -  x(H2S))]2  +  K"x(H2S)  f5  -  0.5(1  -  x(H2S)  +  2K"x(H2S))  —  (A.83) 

and  p(H2S)  =  wm  /  H'(H2S)(1  -  K"). 

The  corresponding  result  for  the  case  x(H2S)  =  0  can  be  derived  as  follows: 

Define  K,  =  (K,KIYH'(C02)  /  Kw)  W  /  y2) 
and  K2  =  (KmKw  /  K,K1YH'(C02))  a 
•••  K,K2  =  KMaV  /  y2 

In  this  notation,  Equations  (A.  15)  and  (A.  16)  become 
p(C02)  =  (1  /  KiK2)  zB  /  (m  -  z  -  B)2 
0  =  m  -  z  -  B  -  z  /  [K2(B  -  z)] 

To  eliminate  z  from  these  equations,  expand  to  get 
z2  -  z(m  +  (1  /  K2))  +  B(m  -  B)  =  0 

•••  z  =  0.5(m  +  (1  /  K2))  ±  (  [0.5(m  +  (1  /  K2))]2  -  B(m  -  B)  f5 
When  m  =  0,  z  =  0  and  B  =  0  (i.e.  no  C02  present), 

0  =  (1  /  2K2)  ±  ([1  /  2K2]2)0'5 

As  K2  >  0,  the  sign  of  the  discriminant  must  be  negative. 

z  =  m/2  +  1/2K2  -  ([B  -  m/2  -  1/2K2]2  +  B/K2)0'5 
Now  p(C02)  =  (1  /  KiK2)  zB  /  (m  -  z  -  B)2 
and  m  -  z  -  B  =  z  /  K2(B  -  z). 

•••  P(C02)  =  (B  /  Ki)-(z  -  B)  /  (z  -  m  +  B) 

Substituting  the  expression  for  z, 

p(C02)  =  (B  /  Ki))  (Q  +  0.5m  -  B)  /  (Q  -  0.5m  +  B)  —  (A.84) 

where  Q  =  (1  /  2K2)  -  ([B  -  0.5m  -  0.5  /  K2]2  +  B  /  K2f5 

But  B  =  mx(C02)  -  H'(C02)p(C02).  These  two  equations  can  be  solved  for  B  and  p(C02). 
Determination  of  K\\ 

When  x(H2S)  =  0,  x(C02)  =  0.5,  Equations  (A.  18)  and  (A.  19)  become 
p(C02)  =  0.5y2  /  (KMa2a2)  (0.5  -  [HC031  /  m)  /  ([HCCEl  /  m)2 
[HCOjI  /  m  =  ([Ki  /  2am)]2  +  K,  /  2am)05  -  K,  /  2am  where  K,  =  K,K1YH'(C02) 

/  KwKm 

Substituting  this  expression  in  the  equation  for  p(C02), 
p(C02)  =  0.5Kwmy2  /  K,K1Y  H'(C02)aa2 

A  K,y  =  0.5my2Kw/p(CO2)aa2K,H'(CO2)  —  (A.85) 

Tertiary  amines  do  not  form  carbamates.  An  equation  similar  to  Equation  (A. 24)  can  be 
obtained  for  moderate  loadings  by  neglecting  [H^,  [OFH,  [C02]  and  [C03  ]  in  the  balance 
equations. 

[HC031  =  mx(C02)  -  [C02]  =  B 
[RR'NH21  =  [HC031  =  B 
[RR'NH]  =  m  -  [RR'NH2+]  =  m  -  B 


. 
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Multiplying  Equations  (A. 5),  (A.  11),  (A.  13)  and  dividing  by  Equation  (A. 6), 

KiK,YH'(C02)  /  Kw  =  (y2/aa2)  B2/[(m  -  B)p(C02)] 

As  [C02]  «  [HC03~],  B  =  [HCCEl  as  mx(C02). 

•••  Kiy  =  mx2(C02)y2Kw  /  p(C02)(l  -  x(C02))att2K,H'(C02)  —  (A.86) 

Determination  of  Km  for  primary  and  secondary  amines 

For  moderate  loadings,  Equations  (A.  18)  and  (A.  19)  hold  with  x(H2S)  set  equal  to  zero.  These 
can  be  solved  for  Km  by  trial  and  error.  (Assume  KM.  Calculate  Ki,  [HCC^l/m,  p(C02)  and 
check  against  the  experimental  p(C02).  Repeat  till  the  predicted  and  the  experimental  p(C02) 
are  same.) 

Determination  of  K\ c 

Consider  the  system  H2S  /  amine  /  water.  For  moderate  loading.  Equation  (A.  17),  with 
x(C02)  =  0,  gives 

p(H2S)  =  (Kw/K,K,cH'(H2S))-  (y2/aa)-mx2(H2S)  /  (1  -  x(H2S)) 

Kic  =  (Kw/p(H2S)K,H'(H2S))-  (72/aa)  mx2(H2S)  /  (1  -  x(H2S)) 
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APPENDIX  B 


Brown’s  Method 

In  this  Appendix,  Brown’s  method  of  solving  a  system  of  nonlinear  equations  is 
described.  The  method  has  been  discussed  in  detail  by  Brown  (1969),  Brown  and  Dennis 
(1971)  and  Brown  (1973).  Brown’s  method  is  based  on  a  partial  pivoting  technique  similar  to 
Gaussian  elimination. 

Consider  the  following  system  of  algebraic  equations: 
fi(x)  =  f i( x i ,  x2,  .  .  .  ,  xN)  =  0 
f:(x)  =  f2(xi,  x2,  .  .  .  ,  xN)  =  0 

fN(x)  =  f nt( x i ,  x2,  .  .  .  ,  xN)  =  0 
or,  in  vector  form,  f(x)  =  0. 

In  Newton’s  method,  the  functions  are  expanded  simultaneously  about  a  point  xn 
assumed  to  be  close  to  the  solution  x*. 

f(x)  =  f(xn)  +  J(xn)-(x  -  x")  +  Higher  order  terms 
J  is  the  Jacobian  matrix  defined  by  Jy  =  3f,/3xj.  When  x  is  close  enough  to  xn,  the  higher 
order  terms  are  negligible.  If  it  is  assumed  that  x  is  a  close  approximation  to  x*,  so  that  f(x) 
=  0, 

0  =  f(xn)  +  J(xn)  -(x  -  xn) 

Solving  for  x, 

xn+1  =  x  =  xn  -  J~‘(xn)-f(xn) 

Thus  starting  with  a  guess  x°,  one  can  solve  iteratively  for  n  =  0,  1,  2,  ...  .  When  the  partial 
derivatives  3f/3xj  are  estimated  by  the  first  difference  approximation,  the  method  is  called  the 
discrete  Newton’s  method.  In  this  method, 

3f,/ 3 Xj  =  [f,(xn  +  hnej)  -  f(xn)/hn 

where  e,  is  the  j-th  unit  vector  and  hn  is  chosen  so  that  hn  ||  f(xn)  ||. 

Brown’s  method  handles  the  functions  one  at  a  time  and  uses  the  information 
obtained  from  working  with  fi  while  working  with  f2,  etc..  This  successive  substitution  nature 
of  Brown’s  method  is  similar  to  Gaussian  elimination  for  linear  systems.  Numerical  stability  of 
this  method  is  improved  by  partial  pivoting  as  described  below. 

Brown’s  method  requires  fewer  function  evaluations,  (N2+3N)/2,  as  compared  to 
Newton’s  method  which  requires  (N2+N)  evaluations  in  each  step.  The  method  is  at  least 
second  order  convergent.  This  means  that  the  number  of  significant  digits  approximately 
doubles  in  each  iteration.  The  method  consists  of  the  following  steps: 

Step  1.  Expand  the  first  function  fj  in  a  Taylor  series  expansion  around  xn.  Ignoring  the 
higher  order  terms, 

f,(x)  =  f,(xn)  +  X  fi.xO).h(xn)-(xJ  -  x/) 
where 

fl.x(j).h(Xn)  =  [f,(xn  +  h  Cj)  -  fi(xn)]  /  h 


. 
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Define 

f."  =  f,(xn) 

and 

fl"x(j>-h  =  fl.x(j).h(x") 

If  x  is  close  enough  to  x*,  fi(x)  =  0.  Thus  we  can  solve  for  that  variable,  say  Xn,  whose 
corresponding  derivative  finx(N).h  is  the  largest  in  absolute  value. 

Xn  =  Xn  -  (f l n/fi nx(N).h)  -  2  (fl "x(j).h /f  1  nx(N).h) '( Xj  -  Xj”) 

The  constants  (fi n/ fi nX(N).h)  and  (finxo)h/finx(N).h)  are  stored  for  future  use.  Here  xN  is  a  linear 
function  of  the  N-l  variables  xi,  X2,  .  .  .  ,  xn-i.  Rename  the  left  hand  side  as  Ln(xi,  X2,  .  .  . 

,  Xn-i) 

Xn  =  Ln(Xi,  X2,  .  .  .  ,  Xn— i). 

Also  define  LN"  =  LN(xin,  x2n,  .  .  .  ,  Xn-i"). 

Step  2.  Define  a  function  g2  of  the  N-l  variables  xi,  x2,  .  .  .  ,  xN_i  as  follows. 
g2(Xi,  .  .  .  ,  Xn-i)  =  f2(Xi,  .  .  .  ,  Xn-1,  Ln(Xi,  .  .  .  ,  Xn-i) 

and 

g2  —  l2(Xi  ,  .  .  .  ,  XN— 1  ,  Ln  ) 

Now  expand  g2  in  a  Taylor  series  and  solve  for  that  variable,  say  xN_i,  whose  corresponding 
derivative  g2.X(N— i >.h  is  the  largest  in  magnitude. 

Xn— 1  =  Xn— r  -  (g2n/g2nx(N— l).h)  -  2  (g2nx(j).h/ g2 "x (N— 1  ).h) *( Xj  -  Xj") 

Define  the  right  hand  side  as  Ln-i(xi,  .  .  .  ,  xn -2). 

Step  3.  Define 

gi(Xi,  .  .  .  ,  Xn— 2)  =  f2(Xi,  .  .  .  ,  Xn— 2,  Ln— 1,  Ln) 
where  LN-i  =  LN-i(xi,  x2,  .  .  .  ,  xN->) 
and  Ln  =  Ln(xi,  x2,  .  .  .  ,  Xn— 2,  Ln— 1). 

The  process  in  Steps  1  and  2  is  repeated  and  the  variable,  say  Xn_2,  whose 
corresponding  derivative  g3.x(N-2).h  is  largest  in  magnitude  is  written  as 
Ln— z(xi,  X2,  .  .  .  ,  Xn— 3). 

In  this  way,  each  step  adds  one  linear  equation  to  the  linear  system.  During  (k+l)-st 
step,  it  is  necessary  to  evaluate  gk+i  (—  fk+i)  for  various  arguments.  These  arguments  are 
obtained  by  back-substitution  in  the  system  Ln,  Ln— 1,  .  .  .  ,  LN-k  which  has  built  up. 

Step  N.  Here  we  have  gN  =  fN(xi,  L2,  .  .  .  ,  LN) 

which  is  a  function  of  the  single  variable  Xi.  Expand  gN,  linearize,  and  solve  for  xi. 

Xl  —  X|  -  gN  /  gN  x ( 1  ).h 

The  point  X]  obtained  from  this  expression  is  used  as  the  next  approximation  Xi"+1 
to  the  first  component  Xi*.  The  complete  Lj  system  is  back-solved  to  get  improved 
approximation  to  X2*,  .  .  .  ,  xN*. 

Here 

Xjn+1  =  Lj(Xln+1, . . . ,  Xj_r+l) 

The  details  on  convergence  theorems,  choice  of  h  etc  may  be  found  in  the  references  cited 
above. 


65 


APPENDIX  C 


Computer  program 

The  computer  program  described  in  Chapter  IV  is  available  upon  request  from  the 
Department  of  Chemical  Engineering,  The  University  of  Alberta,  Edmonton,  Alberta,  Canada. 
A  flow  chart  of  this  program  is  presented  on  the  next  page. 


START 


Read  the  temperature,  amine  code 
and  amine  normality 


Calculate  the  equilibrium 
constants,  Henry’s  constants 
and  dielectric  constant  of  H2Q 


Read  the  partial  pressures 
or  the  loadings  of  acid  gases 


Calculate  a  preliminary 
solution  assuming  7 *  =  1 
and  </>,  =  1 

Refine  the  p 
estimate  using  B 

reliminary 
rown’s  method 

Print  the  loadings 
or  the  partial  pressures 


C 


STOP 


3 


Figure  28.  Flow  chart  for  the  computer  program 


